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Abstract 

Branching fractions and CP-violating asymmetries of charmless B s — > PP, VP, VV decays (P and V 
denoting pseudoscalar and vector mesons, respectively) are re-examined in the framework of QCD factor- 
ization (QCDF). We take into account subleading power corrections to the penguin annihilation topology 
and to color-suppressed tree amplitudes that are crucial for resolving the CP puzzles and rate deficit prob- 
lems with penguin-dominated two-body decays and color-suppressed tree-dominated n Q K° and p 7T° modes 
in the B u j sector. Many of the B s — > h\\i2 decays can be related to B c / — > h\h2 ones via t/-spin or SU(3) 
symmetry. Some useful model-independent relations can be derived and tested. Mixing-induced CP asym- 
metries for many of the penguin-dominated decays are predicted to be very small in the standard model. 
They are sensitive to New Physics and offer rich possibilities of new discoveries. Measurements of direct 
CP-violating asymmetries can be used to discriminate QCDF from other competing approaches such as 
pQCD and soft-collinear effective theory. 
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I. INTRODUCTION 



The phenomenology of nonleptonic two-body decays of B mesons offers rich opportunities for our un- 
derstanding of the underlying mechanism for hadronic weak decays and CP violation. In the past decade, 
nearly 100 charmless decays of B u j mesons have been observed at B factories, BaBar and Belle, with a 
statistical significance of at least four standard deviations (for a review, see ill]). The CDF Collaboration has 
made unique contributions to the measurements of charmless hadronic B s decays. Recently, Belle has also 
started to study the weak decays of the B s meson. 

Many of the B s — > hiti2 decays can be related to Bd — > h[h' 2 ones via £/-spin or SU(3) symmetry. Some 
useful model-independent relations can be derived and tested. For example, direct CP asymmetries of 
B s — > K + n~ and B c i — > K~n + are related to each other by £/-spin symmetry. Therefore, the use of flavor 
symmetry will be helpful to control the hadronic uncertainties in B s — > h\ti2 decay amplitudes. 

Analogous to the neutral Bd system, CP violation in B s decays also occurs through the interference of 
decay amplitudes with and without B s — B s mixing. It is known that the mixing-induced CP violation of 
Bd —* J/yK is governed by sin 2/3. Likewise, the decay B s — ► J/yQ is the benchmark in the B s system 
with mixing-induced CP asymmetry characterized by sin2/3 s . In the standard model (SM), the phase /3 4 
is very small, of order 1 degree. Consequently, B s — > J/y<p and several charmless penguin-dominated B s 
decays e.g. B s -» K^°K^° , rj rj ^ , are the ideal places to search for New Physics as CP violation from 
physics beyond the SM can compete or even dominate over the small SM CP phase. Recently, both CDF 
|Ql and DO ^\ have observed 1-2 a deviations from the SM prediction for f5 s . Because of the possibilities 
of new discoveries, the search for New Physics in the B s system will be the main focus of the forthcoming 
experiments at Fermilab, LHCb and Super B factories. 

Theoretically, two-body B, decays have been studied in the framework of g eneralized factorization 
QCD factorization (QCDF) flflflflS], perturbative QCD (pQCD) 0, [ill 13 and soft-collinear effective 
theory (SCET) Ill3l. ll4ll. In this work we will re-examine and update the QCDF predictions. Especially, we 
shall pay attention to the issue of power corrections. From the study of charmless hadronic B u j decays, we 
learned that two subleading 1 jm\, power corrections are needed in QCDF in order to account for the observed 
rates and CP asymmetries. Power corrections to the penguin annihilation topology, corresponding to the so- 
called "scenario S4" in |gj, are crucial for accommodating the branching fractions of penguin-dominated 
B u ,d PP,VP,VV decays on the one hand and direct CP asymmetries of Bd — ► K~n + , Bd — > K*~n + , 
B — > K~p° and Bd — ► 7T + 7T~ on the other hand. Otherwise, the predicted rates will be too small and 
CP-violating asymmetries of above-mentioned modes will be wrong in signs when confronted with exper- 
iment. However, power corrections due to penguin annihilation will bring new CP puzzles for the decays 
B -^K~K { \ K~r\,n~r\, Bd — > K*°r\ and Bd —> K°7i°: Signs of their A C p's are flipped into the wrong ones 



when compared with experiment. It has been shown in H 1 3TI that soft corrections to the color-suppressed tree 
amplitude due to spectator scattering and/or final-state interactions will bring the aforementioned CP asym- 
metries to the right track and accommodate the observed 7i°7t° and p°7t° rates simultaneously. 1 Recently 



It is well known that a large complex electroweak penguin can also solve the B — > Kn CP puzzle with the differ- 
ence of Acp{B~ — > K^u ) and Acp(B° — > K~n + ) (see e.g. |16]). Since the electroweak penguin amplitude /few is 



2 



we have given a detailed study of charmless hadronic B ur j — > PP, VP, VV decays within the framework of 
QCDF incorporating aforementioned power coiTections I17I1 . In this work we shall generalize the study to 
B s decays. So far B s — ► K + 7t~ is the only hadronic decay mode in the B s sector that its direct CP violation 
has been measured 11811 . The resulting CP asymmetry Acp{B s — ► K + n~) = 0.39 ±0.17 differs from zero 
by 2.2a deviations. Just as the decay B c \ — > K~7i + , the predicted CP asymmetry for B s — > K + K~ in the 
heavy quark limit is wrong in sign and too small in magnitude. As we shall see below, we need penguin 
annihilation to get the right sign and magnitude for Acp(B s — ► K + n~). 

This work is organized as follows. We outline the QCDF framework in Sec. 2 and specify various 
input parameters, such as form factors, light-cone distribution amplitudes and the parameters for power 
corrections in Sec. 3. Then B s — > PP, VP, VV decays are analyzed in details in Sees. 4, 5 and 6, respectively. 
Conclusions are given in Sec. 7. 



II. B DECAYS IN QCD FACTORIZATION 



Within the framework of QCD factorization 11911 . the effective Hamiltonian matrix elements are written 
in the form 

(M l M 2 \^ eS \B) = ^ £ X^{M { M 2 \3?J^ P +^ P \B), (2.1) 

where Xp^ = V p bV* q with q = d,s, and the superscript h denotes the helicity of the final-state meson. For 
PP and VP final states, h = 0. £? £ / h,p describes contributions from naive factorization, vertex corrections, 
penguin contractions and spectator scattering expressed in terms of the flavor operators af ,h , while ^ 



contains annihilation topology amplitudes characterized by the annihilation operators bf ,h . Specifically 1 19] 

S?J X = a p l (MiM2)8p U (ub) v _ A ®(qu) v _ A +aP(M l M 2 )8 pu {qb) v _ A ®{au) v _ A 
+ of (M 1 M 2 )Y,iqb) v _ A 8) (q'q') v . A +a p 4 (M l M 2 )Y,(<i'b)y_ A ® (qq') v _ A 
+ a p 5 {M y M 2 )^(qb) v _ A ® (q'q') v+A +a p (M l M 2 ) J £(-2)(q'b) s _ p ® (qq') s+p (2.2) 

+ of {M x M 2 )^{qb) v _ A ® ^e g (q , q')y +A +a p (M 1 M 2 )'£(-2)(q'b) s _ p ®^(qq') s+ p 
+ a p {M Y M 2 )^{qb) v _ A O ^e g (q'q , ) v _ A +a p Q (M 1 M 2 )'£(q'b) v _ A ® \e q (qq') v _ A , 

where {qiq 2 ) v±A = <7i7/i(l ± 7s)^2 and (qiq 2 ) S:kP = ^i (1 ± 75)^2 an d the summation is over q' = u,d,s. The 
symbol (8) indicates that the matrix elements of the operators in 5^ are to be evaluated in the factorized 



essentially real in the standard model, one needs New Physics to produce new strong and weak phases for Pew- In 
principle, it will be difficult to discriminate between large complex color-suppressed tree C and large Pew scenarios 



in the penguin-dominated decays. However, as pointed out in 11711 . the two schemes can lead to very distinct predic- 
tions for tree-dominated decays where Pew C. The observed decay rates of B° — ► 7r 7r°, p°7r° and the CP puzzles 
with 7t~i~i and 7r°7r° can be resolved by a large complex C but not Pew- In the B u c / sector, there are 13 modes in 
which CP asymmetries have been measured with significance above 1.8<T: K~n + , n + n^ rj ,K* Q rj p° ,p ± n zf 



and p+ K ,K*~ % + ,K n° ',k~ 7] ,coK° ',n c 'ft ', p~ k + . We have shown in ll7|] that the QCDF predictions of A C p for 
aforementioned 13 decays are in agreement with experiment except the decay B° — ► a>K°. However, we notice that 
BaBar and Belle measurements of Acp((oK°) are opposite in sign. 
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form. For the decays B — > PP,VP, VV, the relevant factorizable matrix elements are 

= (P 2 |y A( |0)(P 1 |^|B)=/ /p2 (ml-m2 i )F ( f Pl (O, 
X (UP,V) _ (v|7 M | )(P|^|B)=2/ y m BA .Ff p (m2), 
x (Bv,p) _ ( P | 7 |0)(y|/M|B) =2 / pmBPc Af 



(BVuVz) 



= (V 2 \J fl \0}(V 1 \J , ^\B) = -if V2 m 2 



(et-4){m B + m Vl )A* Vi (r4 2 ) 



(e 1 *-p fl )(e 2 *-p B ) 



2A 



W2 



J-ic C*P C* V n a rfi 



2 yBV i ( 



m 



(rriB + m Vl ) 



(2.3) 



where we have followed the conventional definition for form factors [20]. For B — > VP,PV amplitudes, 
we have applied the replacement my£* • p B — > m B p c with p c being the cm. momentum. The longitudinal 
(h = 0) and transverse (h = ±) components of xjf^ 1 '^ are given by 



X 



(BV h Vz) 



v 2 







2m 



[m l B - m\ x - m l Vi ) (m B + m Vl )Af Vl (q 



X 



(BV U V 2 ) 



-ifv 2 m B m V2 



1 + 



.BV, 



(* 2 )T 



2p c 



m B + m Vl 



4m * P " AfV) 
V BV '(<7 2 



(2.4) 



The flavor operators af' A are basically the Wilson coefficients in conjunction with short-distance non- 
factorizable corrections such as vertex corrections and hard spectator interactions. In general, they have the 



expressions 



af' n (MiM 2 ) 



V? (M 2 ) + ^H'/iM^)] + Pf > P (M 2 ), (2.5) 



N c J-' ' N c An 

where i = l,- - ,10, the upper (lower) signs apply when i is odd (even), c; are the Wilson coefficients, 
Cp = — l)/(2N c ) with N c = 3, M 2 is the emitted meson and M\ shares the same spectator quark with the 
B meson. The quantities v! l (M 2 ) account for vertex corrections, Hj 1 (M\M 2 ) for hard spectator interactions 
with a hard gluon exchange between the emitted meson and the spectator quark of the B meson and Pi(M 2 ) 
for penguin contractions. The expression of the quantities Nf(M 2 ) reads 



N?{M 2 




(2.6) 



The weak annihilation contributions to the decay B — > M\M 2 can be described in terms of the building 

(2.7) 



blocks bf ,h and bf£ w 



% £ 4^(M 1 M 2 \£r^p\^ 

V2 , 



V ^ p=u,c i 



p=u,c 

The building blocks have the expressions 

Cf 



i Cf 

3,EW = AT? 

C F 



c 2 A\, 



b. 



Ni 
C F 

f\ 



c 3 A\ + c 5 {A 3 +a{) + N cC6 A{ 



caA\ + C6A' 



c 9 A\+ Cl {A\ +Ai i ) + N c c%A\ 



£>4,ew = [ c wA\ + csA' 2 ] . 



(2-8) 
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Here for simplicity we have omitted the superscripts p and h in above expressions. The subscripts 1,2,3 of 
A l n denote the annihilation amplitudes induced from (V -A)(V — A), (V -A)(V+A) and (S — P)(S + P) 
operators, respectively, and the superscripts i and / refer to gluon emission from the initial and final-state 
quarks, respectively. Following [a] we choose the convention that Mi contains an antiquark from the weak 
vertex and M 2 contains a quark from the weak vertex. 

For the explicit expressions of vertex, hard spectator corrections and annihilation contributions, the 
reader is referred to for details. The decay amplitudes of B s — > PP, VP are given in Appendix 

A of [6] and can be easily generalized to B s — > VV (see Jsfl for explicit expressions of B s — > VV amplitudes). 
In practice, it is more convenient to express the decay amplitudes in terms of the flavor operators af ,p and 
the annihilation operators fif which are related to the coefficients a h i ,p and bf by 

af(MiM 2 ) = a\{M x M 2 ), 
d{(M x M 2 ) = a\(MiM 2 ), 

a\' p (M,M 2 ) = / a 'i P{MllV1 



af (MjM 2 ) - a' 1 / (M X M 2 ) for M X M 2 = PP, VP, 
(M X M 2 ) + a 5 p (MiM 2 ) for M X M 2 = VV, PV, 



< P (M X M 2 ) = | t (MlM2) + f f ° rM ^= PP > PV > (2.9) 

' a h 4 ' p (M\M 2 ) - rf a\ p {M X M 2 ) for M X M 2 = VP , V V, 

a h / (M X M 2 ) - a^ p (M X M 2 ) for M { M 2 = PP, VP, 
K a/ (M X M 2 ) + a 1 / (M X M 2 ) for M X M 2 = VV, PV, 

Kp (mm) I a% p (M 1 M 2 ) + r^a h /(M 1 M 2 ) for M\M 2 = PP, PV, 
W 1 2> ' \ t ^(M l M2)-r^4 P {M 1 M 2 ) for M l M 2 = VP,VV, 



a^ w (M 1 M 2 ) 



and 



fif(M\M 2 ) = lfB l M ^ M l b P (2.10) 

The order of the arguments of a p (MiM 2 ) and j3f (M[M 2 ) is consistent with the order of the arguments of 
xiBMiMi) =a m ^_ c hhal factor r x is given by 

r zW = — 7 — 77 1 77^:, r vW = — — 7 — • (2-11) 

HI. INPUT PARAMETERS 



It is clear from Eq. (12.31 ) that we need the information on decay constants and form factors in order 
to evaluate the factorizable matrix elements of 4-quark operators. Moreover, we also need to know the 
light-cone distribution amplitudes of light hadrons in order to evaluate the nonfactorizable contributions. 

A. Form factors 

There exist one lattice and several model calculations of form factors for B s — > P, V transitions: 
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1). In the pQCD approach, the relevant form factors obtained at q 2 = are 111 111 (for simplicity, form 
factors hereafter are always referred to the ones at q 2 = 0, unless specified otherwise) 



-B S K _ n 9 4 +0.05+0.00 
~~ u ' ZH -0.04-0.01 ' 

tB s K* 



0.21 



+0.04+0.00 
-0.03-0.01' 



V 



n 9^+0.05+0.00 
u - ZJ -0.04-0.01' 



r 

I B S K* 




0.30 



+0.06+0.01 



0.25 



0.05-0.01 ' 

+0.05+0.00 
-0.05-0.01' 



n o n +0.05+0.00 
u - JU -0.05-0.01 ' 



iB.K* 



{ B S 4 



n i £+0.03+0.00 

- u - iD -0.03-0.01' 

ft i Q+0.03+0.00 
u,l: '-0.03-0.01- 



(3.1) 



2). Form factors obtained by QCD sum rules are 



7 b s k 



0.30 



+0.04 
-0.03' 



for the B s —> K transition 12111 and 



V B ^ = 0.311 ± 0.026, Af** = 0.360 ± 0.034, 



V 1 



0.434 ±0.035, 



0.474 ±0.033, 



for B s ->• V transitions [22]. 



A?*** =0.233 ±0.022, 
A, 50 =0.311 ±0.030, 



(3.2) 



(3.3) 



3). Another light-cone sum rule calculation based on heavy quark effective theory gives H23I1 



F*' K = 0.296 ±0.018, 



j?B s T] _ n 9 o 1+0.015 



(3.4) 



and 



\/B s K* _ r, 90^+0. 013 \B S K* _ r, 9^9+0.011 \B S K* _ r, 99-y-t-u.uiu 
v — u - zo;, -0.013' "^0 — u - zzz -0.010' A \ ~ u - zz/ -0.012' 



^=0.339+°^ 



\B S $ 



269+ 014 



-0.017' 



\B>4 



0.271 



+0.010 
-0.012 

+0.014 
-0.014- 



(3.5) 



It is clear that form factors obtained by sum rules are larger than the pQCD ones. 



4). A light cone quark model in conjunction with soft collinear effective theory was constructed in H24TI . 
The predictions are 



Fq' k = 0.290, pB s r\ s = Q 288 ^ 
V B * K * =0.323, A BsK * = 0.279, A 



0.228, 



v B.t = 0.329, aJ^ = 0.279, A B = 0.232 . 



(3.6) 



5). A straightforward application of the covariant light-front quark model of M25H yields M2611 



V BsK * =0.23, 
V B '* = 0.30, 



0.26, 
0.32, 



A BsK * =0.19, 
A? 50 =0.26, 



(3.7) 



all with errors estimated to be ±0.01 . 



6). A recent lattice QCD calculation yields F B sK = 0.23 ± 0.05 ± 0.04 1271] 



For comparison, Beneke and Neubert |6j] used 



^ = 0.31 ±0.05, 



l 



0.29 ±0.05, 



0.34 ±0.05, 



(3.8) 



and 



n B,^7) ( 



f q 







+ F 7 



(3.9) 



while Beneke, Rohrer and Yang [8.] employed 

a BsK' 



0.33 ±0.05, 



0.38 



-o.io 



q — <J.^ -J-XJ.xju, ziq — xj.^u _ _02- (3.10) 

Note that it is most convenient to express the form factors for B—*r\0 transitions in terms of the flavor 
states qq = (uu + dd) / y/2, ss and cc labeled by the T] q , t] s and T]°, respectively. Neglecting the small mixing 
with T7 C °, we have 

F B S 7] = _pB s 7] s sm Q F B S 7)' = F B S 7) S CQsd ^ 

where 6 is the r] q — T] s mixing angle defined by 

| T7> = cos0|T7 9 )-sin0|T7 i .), 



(3.11) 



|tj') = sin0|7] ? ) +cosQ\r\ s 



(3.12) 



with 6 = (39.3 ± 1.0)° in the Feldmann-Kroll-Stech mixing scheme |28|]. 

From the above discussions we see that the form factor F Q 1 at q = ranges from 0.23 to 0.31 . In the 



QCDF approach, if F Q * (0) = 0.31 is employed, we find that the predicted branching fractions S3{B S — > 
K + k~) w 9.1 x 10~ 6 and £8{B S — > K + K~) « 34 x 10 -6 will be far above the experimental measurements of 



(5.0 ± 1.1) x 10" 6 |29j] and (25.7 ±3.6) x 10~ 6 13(1 Bill, respectively. Hence we shall use F^(0) = 0.24 
obtained by the lattice calculation. Note that a % 2 analysis by one of us (C.K.C.) with the available data of 



B s -» PP also yields Fq sK (0) = 0.240^^^ [QJ. For other form factors, we shall use F^ s ' h (0) = 0.28 and 
B s — ► V transition form factors given by Eq. (13.71 ) with some modifications on B s — ► K* ones (see Table U) . 



B. Decay constants 



Decay constants of various vector mesons defined by 

(V(p,e)\q 2 Yfiqi\0) = -j/ymye*, 
(V(p,e)\q 2 Opvqi\0) = -ri(s; P v -e* v pn, (3.13) 
are listed in Table J] They are taken from H34n . For pseudoscalar mesons, we use f n = 132 MeV and 



/k = 160 MeV. Decay constants f^^, f s ^ } and f^ {l) defined by 



(0|?7 M 75?ln 



On 



i-^f n( > } q^ (^r l iY5s\r] ( ' ) )=if r , ( ' ) q^ {0|cy M 7 5 c|7](')) = if^q^ 

77 O ' '77O' 



V2^ ( 

are also needed in calculations. For the decay constants f q , and f s ,,,, we shall use the values 



4 = 107MeV, /*=-112MeV, jfj = 89MeV, = 137 MeV 



obtained in 12811 . As for f^, a straightforward perturbative calculation gives 13511 



2 iA 

111 (l\ J (l\ 

7)' ) J 7/( > 

12^2 7f 



(3.14) 



(3.15) 



(3.16) 
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TABLE I: Input parameters. The values of the scale dependent quantities fy (ju) and a\j (m) we given for 



jU = 1 GeV. The values of Gegenbauer moments are taken from 13611 and Wolfenstein parameters from 13711 . 



Light vector mesons 


V 


/y(MeV) 


#(MeV) 


a] 7 




_L.V 
a l 


±.v 

a 2 


P 
ft) 

K* 


216±3 
187 ±5 
215±5 
220 ±5 


165 ±9 
151 ±9 
186±9 
185 ±10 






0.03 ±0.02 


0.15 ±0.07 
0.15 ±0.07 
0.18±0.08 
0.11 ±0.09 






0.04 ±0.03 


0.14 ±0.06 
0.14 ±0.06 
0.14 ±0.07 
0.10 ±0.08 


Light pseudoscalar mesons 


< 


a£ 


af 


af 





0.25±0.15 


0.06 ±0.03 


0.25 ±0.15 


B mesons 


B 


m B (GeV) 


T B (ps) 


/g(MeV) 


Afi(MeV) 


B u 


5.279 


1.638 


210 ±20 


300 ±100 


B d 


5.279 


1.525 


210 ±20 


300 ±100 


B s 


5.366 


1.472 


230 ±20 


300 ±100 


Form factors at q 2 = 


F BsK (0) 


A^*(0) 


Af s **(0) 


A^(O) 


v BsK \o) 


0.24 


0.30 ±0.01 


0.24 ±0.01 


0.22 ±0.01 


0.28 ±0.01 




A B /(0) 


A? S0 (O) 


A*'*(0) 




0.28 


0.32±0.01 


0.26 ±0.01 


0.23 ±0.01 


0.30 ±0.01 


Quark masses 


mj,(mb) /GeV 


m c{ m b) /GeV 


pole / pole 

ml /ml 


m s (2.l GeV) /GeV 


4.2 


0.91 


0.3 


0.095 ±0.020 


Wolfenstein parameters 


A 


A 


P 




7 


0.8116 


0.2252 


0.139 


0.341 


(67.8^|)° 



8 



C. LCDAs 



We next specify the light-cone distribution amplitudes (LCDAs) for pseudoscalar and vector mesons. 
The general expressions of twist-2 LCDAs are 



<t> P (x,lJ,) = 6x(l -x) 



&Ux,n) = 6x(l-x) 



&[(x,/j.) = 6x(l -x) 



l + £a P Md /2 (2x-l) 

n=i 

CO 

\ + ^a v n (^)d /2 (2x-\) 

n=l 

OO 

l + ^at V (^)d /2 (2x-l 



n=\ 



(3.17) 



and twist-3 ones 



%(x) = 1, ® a (x) = 6x(l -x) 



<D v (x,M)=3 



2x-l + ^ay^)P n+l (2x-l) 



(3.18) 



where C n (x) and P n (x) are the Gegenbauer and Legendre polynomials, respectively. When three-particle 
amplitudes are neglected, the twist-3 0> v (x) can be expressed in terms of <I>^ 



u 



du 



1 



du. 



The normalization of LCDAs is 



I dx&y{x) = \, ( dx® v (x)=0. 
Jo Jo 



(3.19) 



(3.20) 



Note that the Gegenbauer moments a • displayed in TableUtaken from yg] are for the mesons containing 
a strange quark. 

The integral of the B meson wave function is parameterized as il9ll 

1 d P ^Bf \ m B 



(3.21) 

lo 1 — p As 

where 1 — p is the momentum fraction carried by the light spectator quark in the B meson. We shall use 
X B = 300 ± 100 MeV. 

For the running quark masses we shall use j38l.l39ll 

m h {m b ) =4.2GeV, m fo (2.1GeV) = 4.94GeV, m,,(lGeV) = 6.34GeV, 

m c (m h ) = 0.91 GeV. m c (2.1GeV) = 1.06GeV, m c (lGeV) = 1.32GeV, 

m,(2.1GeV) =95MeV, m,(lGeV) = 118MeV, 

m d (2.1GeV) =5.0MeV, m„(2.1GeV) = 2.2MeV. (3.22) 



Note that the charm quark masses here are smaller than the one m c {m^) = 1.3 ±0.2 GeV adopted in |g,|9|] and 
consistent with the high precision mass determination from lattice QCD Il40ll : m c (3GeV) = 0.986 ±0.010 
GeV and m c (m c ) = 1.267 ±0.009 GeV (see also j4lll ) Among the quarks, the strange quark gives the major 
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theoretical uncertainty to the decay amplitude. Hence, we will only consider the uncertainty in the strange 
quark mass given by ra s (2.1GeV) = 95 ±20 MeV. Notice that for the one-loop penguin contribution, the 
relevant quark mass is the pole mass rather than the current one Il42ll . Since the penguin loop correction is 
governed by the ratio of the pole masses squared 57 = (m? ole /m p b ° le ) 2 and since the pole mass is meaningful 
only for heavy quarks, we only need to consider the ratio of c and b quark pole masses given by * c « (0.3) 2 . 



D. Penguin annihilation 



In the QCDF approach, the hadronic B decay amplitude receives contributions from tree, penguin, elec- 
troweak penguin and weak annihilation topologies. In the absence of 1 / nib power corrections except for 
the chiral enhanced penguin contributions, the leading QCDF predictions encounter three major difficulties: 
(i) the predicted branching fractions for penguin-dominated B — ► PP, VP, VV decays are systematically be- 
low the measurements, (ii) direct CP-violating asymmetries for B d — > K~n + , Bj — > K*~n + , B — > K~p°, 
Bd — ► 7T + 7T~ and B s — > K + k~ have signs in disagreement with experiment, and (iii) the predicted longitu- 
dinal polarization fractions in penguin-dominated B — ► VV decays are usually too large and do not agree 
with the data. This implies the necessity of introducing 1 /nib power corrections. Unfortunately, there are 
many possible 1 /nib power suppressed effects and they are generally nonperturbative in nature and hence 
not calculable by the perturbative method. 

Power corrections in QCDF always involve troublesome endpoint divergences. For example, the annihi- 
lation amplitude has endpoint divergences even at twist-2 level and the hard spectator scattering diagram at 
twist-3 order is power suppressed and posses soft and collinear divergences arising from the soft spectator 
quark. Since the treatment of endpoint divergences is model dependent, subleading power corrections gener- 
ally can be studied only in a phenomenological way. We shall follow il9ll to model the endpoint divergence 



X = J dx/x in the annihilation and hard spectator scattering diagrams as 



X 



H 



In 



(£)< 



l+p H e" 



(3.23) 



with A/ 7 being a typical scale of order 500 MeV, and Pa,h, <!>a,h being the unknown real parameters. 

A fit to the data of B u j — > PP, VP, PV and VV decays yields the values of Pa and A shown in Table JI] 
Basically, it is very similar to the so-called "S4 scenario" presented in |0]. The fitted Pa and (f>A for B — > VV 
decays are taken from [43]. Since the penguin annihilation effects are different for B — ► VP and B — > PV 
decays, 

r 2~ 



6na s 



3 1 



-VP 



■4 + - 



Jtr x r x 



2X, 



VP 



6na s 
■ 6na, 



rVP\2 



2XY + A 



+ r p J(xD 2 -2Xr + 



n 1 



3r v x (ixr-m 



V VP 
A A 



(3.24) 
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for M\M 2 = VP (the definition for the parameters r£ and rt can be found in Eq. ( 12.111 ) below) and 



■ 6na s 



X 



PV 



(3.25) 



-3^(2XT - l)(2-XD + r\ (2(XD 2 - 

for M\M 2 = PV, the parameters X^ p and X PV are not necessarily the same. Indeed, a fit to the B — > VP,PV 
decays yields pj p « 1.07, <£[ p « -70° and p£ y ~ 0.87, ^ v « -30° (see Table©. For B s . -> PP, VP, 
decays, we shall assume that their default values are similar to that in B u l \ decays as shown in Table JT] For 
the estimate of theoretical uncertainties, we shall assign an error of ±0.1 to Pa and ±20° to <p A . 



TABLE II: The parameters p A and §a for penguin annihilation. 



Modes Pa 0a 


Modes Pa <Pa 


B^PP 1.10-50° 
B^VP 1.07 -70° 
B^PV 0.87 -30° 
B -> K*p 0.78 -43° 
B^K*<p 0.65 -53° 


B S ^PP 1.00 -55° 
B s -> VP 0.90 -65° 
B s -» PV 0.85 -30° 
B s VV 0.70 -55° 



E. Power corrections to a% 



As pointed out in 11511 . while the discrepancies between experiment and theory in the heavy quark limit 
for the rates of penguin-dominated two-body decays of B mesons and direct CP asymmetries of B c i — > 
K~n + , B — ► K~p° and — ► 7T + 7r~ are resolved by introducing power corrections coming from penguin 
annihilation, the signs of direct CP- violating effects in B — > K~n°,B~ — ► K r\ and B° — > 7i°7t° are flipped 
to the wrong ones when confronted with experiment. These new B-CP puzzles in QCDF can be explained 
by the subleading power collections to the color-suppressed tree amplitudes due to spectator interactions 
and/or final-state interactions that not only reproduce correct signs for aforementioned CP asymmetries but 
also accommodate the observed Bj — ► n°7i° and p Q n° rates simultaneously. 

Following [15j], power corrections to the color-suppressed topology are parametrized as 

a 2 ^a 2 {l+Pce i * c ), (3.26) 
with the unknown parameters p c and (pc to be inferred from experiment. We shall use <p c ps —70° and 
Pc ~ 1.3 , 0.8 , for B — > PP,VP, VV decays Ill5l. ll7ll. respectively. This pattern that soft power coiTections 
to a 2 are large for PP modes, moderate for VP ones and very small for VV cases is consistent with the 
observation made in 14411 that soft power correction dominance is much larger for PP than VP and VV final 



states. It has been argued that this has to do with the special nature of the pion which is a qq bound state on 
the one hand and a nearly massless Nambu-Goldstone boson on the other hand 

Hi- 



ll 



IV. B s PP DECAYS 

Before proceeding to the numerical results of QCDF calculations, we discuss some model-independent 
flavor symmetry relations in which many of B s — > PP decays can be related to B d — > PP ones by either 
U -spin or SU(3) symmetry. Hence these relations can be used to cross-check the dynamical calculations. 



A. U -spin symmetry 



In the limit of U -spin symmetry, some of B s decays can be related to B d ones. For example, 

A(B s ^K + 7i-) = V: b V ud (K + K-\O u d \B s )+V; b V cd (K + 7l-\0 L d \B s ), 

A(B d K~n + ) = V: b V l<s (K-x + \0't\B d )+V L * h V cs (K-x + \O l ;\B d ), (4.1) 

where the 4-quark operator O s is for the b — ► qqs transition and O d for the b — ► qqd transition. The assump- 
tion of U -spin symmetry implies that under d <-> s transitions, 



(K + 7C-\O u d \B s ) = {K~ K + \O u s \B d ), 
Using the relation 

lm(V: b V ud V cb V: d ) = 
it is straightforward to show that iS 

\A(B S - K+n-)\ 2 - \A(B S - K-n + )\ 2 
and, consequently, 

A CP (B S ^ + 7ZT) = -Acp(Brf - 



-Im(V: b V us V cb V c 
\A(B d - 

3g(B d ^K-7i+) i(B s 



K~iz 



@(Bs 



(4.2) 
(4.3) 
(4.4) 

(4.5) 



From the cuiTent world averages, A CP (B a 



10~ 6 |48|] and the CDF measurement SB(B S 



K~K + 



-0.0981851?, 



= (19.4 ±0.6) x 



(5.0± 1.1) x 10~ 6 1291], it follows that the predic- 



tion Acp(B s — > K + n ) ~ 0.37 under £/-spin symmetry is in good agreement with the experimental result 
0.39 ±0.15 ±0.08 obtained by CDF fea]. Besides A CP (B S -> ^+71"), CDF has also measured direct CP vi- 
olation in the decay B d — > A -- 7T + and obtained 113011 



-0.83±0.41±0.12, 



(4.6) 



which is equal to — 1 under U -spin symmetry. Obviously, the experimental measurement is still limited by 
statistics. 

By the same token, we also have the following U -spin relations 



A C p(B s - 
A CP (B S 
A C p{B s 
A C p(B s - 



K + K~) 



K Q K°) 



K°n°) 



-A C p(B d 
-A C p{B d 
-A C p(B d 
-A C p(B d 



K + K ) 



K Q K° 



K Q 7l°) 



K K 



<%(B d ^%+%-) x(B s ) 



®{B S - 


^K+K-) %(B d ) 


m{B d - 


+ K°K°) r(B s ) 


m[B s - 


>K°K°) T(B d y 


m(B d - 


>K°7i°) T(B S ) 


®(B S - 


>K n°) x{B d y 


. @{Bd - 


-+K+K-) X(B S 



3$(B S 



r(B d ) ■ 



(4.7) 
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Unlike the first £/-spin symmetry relation ( 14.51 ), the above relations cannot be tested by the present available 
data. Nevertheless, they can be checked by our dynamical calculations as shown in Sec.IV.C.5. 



B. SU(3) symmetry 



There are some cases where two-body decays of B d and B s can be related to each other in the limit of 
SU(3) symmetry provided that some of the annihilation effects can be neglected. Let us consider the decay 
amplitudes of the following three pairs in QCDF (60 : 

. k+k-) = £ v; b V pd A Kn (8 pu ai + < + < - « p - - « p 



A(B S - 
A(B d 



p=u,c 



*4,EW + ft 2^ 3 ' EW ^' 



(4.8) 



1 



p=u,c 



and 



A{B S — > K + K 



A{B d ^K--n + ) 



£ VpbVpdAwidpuCCx + < + < EW + Spufa + j8 3 p - -/3 3 p 

;w) 



V + 2^4.Ew)' 



= I v; b v ps 

p=u,c 

+ B RK (8 pu b p + b p 4 +b p 



2^3,ew 2^ EW -' 



1 



I v; b V ps A*K(8 pu ai + < + <ew + H ~ 2^ew), 



(4.9) 



and 



A(B S — > K°K°) 



p=u,c 



a 



1 



-a 



'4 2 4,EW t 



3.EW 



2^4°Ew) 



+ B KK (b p 4 --b 



2"4,EW^ 



p=u,c 



-a p 

2 U 4,EW 



+ « + /3f + j3 3 p EW ), 



(4.10) 



with Ai H k 2 = x( Bshl ' h2 \ where the expressions of the flavor operators a,- in terms of a,- and the annihilation 
operators j8, in terms of bj are shown in Eq. ( 12.9I ). Roughly speaking, Gti is due to the tree topology, 0:4 
comes from the QCD penguin operators O4 and 0(,, oc p EW receives contributions from the electroweak op- 
erators 08 and f?io- From the study of hadronic B ud decays we learn that annihilation effects are negligible 
in tree-dominated modes and dominated by the /?3 term in penguin-dominated decays. Hence, under the 
approximation of negligible annihilation contributions to tree-dominated decays and keeping only the dom- 
inant penguin annihilations in penguin-dominated decays, SU(3) symmetry (or £/-spin symmetry acting on 
the spectator quark of the B meson) implies 11461. 14911 

A(B S K + 7i~) ^A(B d ^ k + k~), A(B S K + K~) mA(B d -^K~n + ), 



A(B X 



K°K°) 



A(B' 



K°7t- 



(4.11) 



As will be discussed later, it turns out that among the relations 

Sg{B s ^K + %-) » m{B d -» 7T+7ZT), A CP (B S - 
3S{B S -» K + K~) » 3${B d -^K-7t + ), A CP (B S - 



-A CP {B d - 
-A CP (B d 
A CP (B- 



>K + K ), 



(4.12) 



the first three ones are experimentally fairly satisfied. 
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TABLE III: CP-averaged branching fractions (in units of 10~ 6 ) of B s — > PP decays obtained in various 
approaches. In the QCDF calculations, the parameters and 0a are taken from Table HT1 pc = 1.3 and 
0c = —70°. Sources of theoretical uncertainties are discussed in the text. The pQCD predictions to LO 
and (partial) NLO are taken from jllll and 111 211 . respectively. For the decays involving an 77 and/or 77', two 
different sets of SCET results are quoted from 11311 . 



Modes 



Class QCDF (this work) pQCD (LO) pQCD (NLO) 



SCET 



Expt. [18, 29] 



t- 




T 


r n+0.4+0.4 
J ' J -0.8-0.5 


7 6+ 33 
'■°-2.5 


6.3; 


^2.6 
-1.9 


4.9±1.2±1.3±0.3 


5.0±1.1 


t- 




C 


1 7+2.5+1.2 
i - ' -0.8-0.5 


16+ 012 

u - lo -0.07 


0.25 


f0.10 
-0.07 


0.76±0.26±0.27±0.17 




t- 




C 


n 7 c + l. 10+0.51 
u - /J -0.35-0.22 


1 1 + 008 


0.08 


f0.03 
-0.02 


0.80±0.48±0.29±0.18 





5? -A 



On' 



B s ->K u ri 



C 



On' 



B s -> 7t u rj 



B s — > K + K 



5? -► ttV 



0^0 



ann 
ann 



2.8 



+2.5+1.1 

-1.0-0.8 



t- 


^K+K- 


p 


25.2+ 


t- 




p 


26.1 + 


t- 


^7777 


p 


10.9j 


t- 


^7777' 


p 


41.2+ 


t- 


-> 77'77' 


p 


47.9+ 



7.2- 9.1 
13.5+12.9 
- 8.1- 9.4 
6.3+5.7 
4.0-4.2 



+27.3+17.8 
-12.9-13.1 



41.6+20.9 
17.1-15.3 



0.05 



+0.03+0.02 
-0.01-0.01 



0.04 



-0.01+0.01 

-0.00-0.00 



0.26 



-0.00+0.10 

-0.00-0.09 
n io+0.0+0.05 
u ' iJ -0.0-0.05 



0.72 



+0.36 
-0.24 



13.6; 

15.6 
8.0 



-8.6 
-5.2 
-9.7 
-6.0 
-5.4 
-3.1 



21.0+Vi 



14.0 



+7.0 
-4.1 



0.05 



+0.02 
-0.02 



0.11 



0.57 



-0.05 
-0.03 



-0.18 



' -0.16 

O 78+0-09 
U ' 28 -0.08 



1.87 



-0.45 
-0.56 



15.6+|> 
18.0+^ 



10.0 



+3.4 
-2.6 



34.9 



-11.6 
- 9.5 



25 ?+ 8 - 3 
zo - z -6.5 



0.03 



+0.01 
-0.01 



0.08 



0.57 
0.29 



-0.03 
-0.02 



-0.24 
-0.22 
-0.12 
-0.12 



0.59±0.34±0.24±0.15 
4.5 ±1.5 ±0.4 ±0.5 
3.9 ±1.3 ±0.5 ±0.4 
18.2 ±6.7 ±1.1 ±0.5 
17.7±6.6±0.5±0.6 
7.1±6.4±0.2±0.8 
6.4±6.3±0.1±0.7 
24.0±13.6±1.4±2.7 
23.8±13.2±1.6±2.9 
44.3 ±19.7 ±2.3 ±17.1 
49.4±20.6±8.4±16.2 
0.014 ± 0.004 ± 0.005 ± 0.004 
0.016 ± 0.0007 ± 0.005 ± 0.006 
0.006 ±0.003 ±0.002+°;°^ 



0.038 ±0.013 ±0.016 



-0.260 
-0.036 



25.7 ±3.6" 



< 1.2 



^This is the average of the CDF and Belle measurements, (24.4 ± 1.4± 3.5) x 10~ 6 |30|] and (38+^±7) x 10~ 6 



13 ill , respectively. The old CDF result on B s — > K + K can be found in 1321. 



C. Numerical results and comparison with other approaches 



We list in Tables [TTT] and [TV] the branching fractions and CP asymmetries of B s — > PP decays evaluated 
in the frameworks of QCD factorization (this work), pQCD to the lowest order (LO) jllll and to the next- 
to-leading order (NLO) il2l1 and soft-collinear effective theory (SCET) 11 3il . For the decays involving an 
77 and/or 77', two different sets of SCET results are quoted from 11311 . corresponding to two distinct SCET 
parameters regarding to the strong phases of the gluonic charming penguin. The expression for the decay 
amplitudes of B s — > PP and VP decays in the QCDF approach can be found in the Appendix of JfJ. 

The theoretical errors in QCDF calculations correspond to the uncertainties due to the variation of (i) the 
Gegenbauer moments, the decay constants, (ii) the heavy-to-light form factors and the strange quark mass, 
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TABLE IV: Same as Table Hill except for the direct CP asymmetries (in %) in B s — ► PP decays. 

Modes Class QCDF (this work) pQCD (LO) pQCD (NLO) SCET Expt [18] 

^^-la-l'l 25.8+^ 20±17±19±5 39±15±8 

36.3+JT^+Jjj 59.4+! 7 2 9 5 88.0±|J -58 ±39 ±39 ±13 

33 - 4 -l:8-2L6 56A -93 96J -2 6 5 -56±46±14±6 

61±59±12±8 

-49.3+f;g+i|;° 19.9 ^ -35.4+^ -14±7±16±2 

37±8±14±4 

- 7 - 7 -L2^:? -23.315;° ~^.6t\i -6±5±6±2 

o.4oi8:ol±o:o4 0A±0A < 10 

-5.0+£j±|;f -0- 6 la5 °- 6 ^ao 7.9±4.9±2.7±1.5 

-1.1±5.0±3.9±1.0 
-0.6+^+^ -L3la2 -°-2lai 0.04 ±0.14 ±0.39 ±0.43 

2.7±0.9±0.8±7.6 

3 - 2 -a6-l:2 L9 -o:5 L4 -o:2 0.9 ±0.4 ±0.6 ±1.9 

-3.7±1.0±1.2±5.6 







T 






c 






c 




->K°ri' 


c 


t- 


-> K + K 


p 


t- 




p 


t- 


^7777 


p 


t- 




p 


t- 




p 







Pew 


n 6 ,+ 1.6+ 1.8 


-0 4+ 03 


40.4+4-0 


t- 




Pew 


a ? q+2.3+31.0 
HZ '^-8. 1-40.9 


20.6+^ 


52.5+Jj 


t- 


■+ tt+tt - 


ann 





i-2+!;I 


0.2+?;? 


t- 


ttV 


ann 





1.2+M 


2+ 01 



and (iii) the wave function of the B meson characterized by the parameter Ag, the power corrections due 
to weak annihilation and hard spectator interactions described by the parameters Pa,h> fa,H, respectively. 
To obtain the errors shown in Tables ITTTllXIVI we first scan randomly the points in the allowed ranges of 
the above nine parameters and then add errors in quadrature. As noted in passing, we assign an error ±0.1 
and ±20° to the default values of and fa, respectively, while p H and fa lie in the ranges < p# < 1 
and < fa < T.%. Specifically, the second error in the table is referred to the uncertainties caused by the 
variation of Pa.h an d fa.H> where all other uncertainties are lumped into the first error. Power corrections 
beyond the heavy quark limit generally give the major theoretical uncertainties. For theoretical uncertainties 
in pQCD and SCET approaches, the reader is referred to the references cited in the table captions. 

1. B s ^K + K-,K°n\K {) r]^ 

As mentioned before, in this work we shall use the form factor Fq ,k (0) = 0.24 obtained by both lattice 
and pQCD calculations. If a larger B s to K transition form factor, say, F^ sK (0) = 0.31, is employed, the 



15 



predicted 3%{B S — > 7^+ 7T ) and ^(B s — > K + K ) will be far above the experimental results. 2 For Fq sK (0) = 
0.24, the calculated ^(fi , — » ^ + 7T~) = (5.3^Qg^Qg) x 10~ 6 is in good agreement with the measurement 



(5.0 ±0.7 ±0.8) x 10 -6 |22]. Notice that although the same value of F s was used in the leading order 



pQCD calculation, a larger branching fraction of order 7.6 x 10 6 was obtained (see Table HIlT). 



Z = -r °bk V ?" 2 ( 4 -!3) 



A recent detailed analysis in 1)5011 indicates that SU(3) and factorization only remain approximately valid 
if the branching fraction of B s — ► K + n~ exceeds its current value of (5.0 ± 1.1) x 10 -6 by at least 50% or if 
the parameter t, denned by 

h F$ n {rn\) m\ -ml 
fx F^ K {ml) m\-m\ 

is more than about 1.2. The analysis goes as follows. Writing the amplitudes A{B~ — ► K°7i~) = V cs V* b P 
and A(B~d — ► 7C~7r + ) = V us V* h Te lS + V cs V* b P, the measured 5~ — > rate sets a constraint on the penguin 
topology P. Since V„fc = |V„fc|e~ !r , the measurement of B c \ — > K~n + will put a constraint on T as a function 
of the unitarity angle y. Under £/-spin symmetry, the amplitude A(B S — > K + 7i~) = V u dV* b T'e iS ' + V c dV* b P' 
can be related to the B rf — > K~K + one by the relations: T = T,P' = P and 5' = 8. The data of B s -> #+7T~ 



will be helpful for pinning down the ratio of P/T. The analysis of 15011 shows that for the value of / to 
be consistent with other determinations and for the strong phases 8 and 8' not different much from each 
other, then either 38(B S — > ^ + 7r~) is at least 50% larger than the current measured value or the parameter t, 
is larger than 1.2 . Our results of § = 1.24 and 38(B S — > K + n~) « 5.3 x 10~ 6 are thus consistent with the 



analysis of 15011 . 

It is known that the predicted direct CP violation for Bj — ► ^"TT" 1 " and 5^ — ► K + 7i~ modes in naive QCDF 
is wrong in sign when compared with experiment (see the predictions in |0]). This discrepancy together with 
the rate deficit problem for penguin-dominated decays can be resolved by introducing power corrections 
coming from penguin annihilation, corresponding to the "S4 scenario" of [6]. Using the values given in 
Table [TT] for the parameters p A and (j) A , we obtain A C p{B d — > K~7t + ) = -(7.41} 5+43)% and A C p(B, 



K + n~) = (20.7+|^l|;|)%, to be compared with the data -0.098+^} 2 I48J] and 0.39 ±0.15 ±0.08 \M, 
respectively. 

The inclusion of soft corrections to the color-suppressed tree topology has two effects: First, it will 
enhance the rates of B s — > K°K°,K°ri by a factor of about 2.5 and B s — > K°r\' slightly. Second, it will flip the 
sign of CP-violating asymmetries of the former two modes. For example, 38(B S — ► K°7i°) is enhanced from 
0.7 x 10~ 6 to 1.7 x 10~ 6 , while A CP (B S -> K°n°) is changed from -0.214 to the order of 0.363 (see Tables 
UTTl and [TVT >. Note that pQCD predictions of branching fractions for the color-suppressed tree-dominated 
decays B s — > K K a ,K°r\^ are much smaller than QCDF and SCET. Nevertheless, pQCD results of A C p's 
for the above three modes agree in signs with QCDF 

We see from Table HVl that SCET predicts a negative sign for Acp(B s — > K°7l°), contrary to QCDF and 
pQCD. This deserves a special discussion. The negative sign of Acp(B s — > K 7T°) has to do with the fact 



that SCET predicts A C p{B d -> K°n°) = (5 ±4±4± 1)% yj]. From the £/-spin symmetry relation (l4~7l) we 
learn that CP asymmetries of B s — > K°n° and B d — > K°7i° are of opposite sign. Although the current world 



2 A larger branching fraction 3§{B S -> K+n-) = (10.2±f;!|) x 1(T 6 was obtained in (6:] within the framework of 
QCDF using the form factor F$ sK (0) = 0.31 ± 0.05 . 
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average A CP (B d -> K Q n ^ 



-0.01 ±0.10 from the BaBar and Belle measurements, -0.13 ±0.13 ±0.03 



I 51H and 0.14 ±0.13 ±0.06 15211 respectively, is consistent with no CP violation, there exist several model - 
independent determinations of this asymmetry: one is the SU(3) relation 15311 



AT(B a 



-AT(B a 



K°K°) 



and the other is the approximate sum rule for CP rate asymmetries 15411 



(4.14) 



(4.15) 



based on isospin symmetry, where Ar(S — ► Kit) = T(B — > Kit) — Y{B — > Kit). This sum rule allows us to 
extract A CP (B d — ► K°n°) in terms of the other three asymmetries in K~ n + ,K~ 7t° ,K°n~ modes that have 
been measured. From the current data of branching fractions and CP asymmetries, the above SU(3) relation 
and CP-asymmetry sum rule lead to -0.073+g$f and -0.15 ±0.04, respectively, for A C />(£ 



An analysis based on the topological quark diagrams yields a similar result —0.08 ~ —0.12 Il55ll . All these 
indicate that direct CP violation should be negative for B d —> K°7t Q and hence positive for B s — ► K°7l Q . 



2. B t ^K + K-,K°K° 



The penguin-dominated decays B s — > K + K ,K°K° have sizable branching fractions of order 25 x 10~ 6 
in QCDF. The corresponding pQCD and SCET predictions are slightly smaller (Table HID). 3 From Eqs. ( 14.91 ) 
and (14.101 ) we see that K + K and K°K° modes differ mainly in the tree contribution (Xi and the annihilation 
term /3i induced by the operator 0\ , both existing in the former but not in the latter. Since these contributions 
are CKM suppressed relative to the penguin terms, the above two modes should have similar rates but 
rather distinct CP asymmetries. Due to the absence of interference between tree and penguin amplitudes, 
CP asymmetry is very small in B s — ► K°K°, less than 1%. Using the world average of Acp{B d — > n + ii~) = 
0.38 ±0.06, .%(B d -> K+K-) = (5. 16 ±0.22) x 10~ 6 fl and SS{B S -»• K+R-) = (25.7 ±3.6) x lO" 6 
|48D, we find from the first £/-spin relation in Eq. d4"77T ) that A C p(B s ->• K + K~) ps -0.077 in the ?7-spin 
limit, which is in excellent agreement with the QCDF prediction. It is very important to measure the direct 
CP asymmetry for this mode. 

In the pQCD approach, direct CP violation of B s — > K°K° vanishes to the lower order as there is only one 
type of CKM matrix element in its decay amplitude, say V,bV t * jllll . To the NLO, penguin loop corrections 
allow other CKM matrix elements enter into the decay amplitude and induce CP asymmetry I12I1 . It turns 
out that the predicted Acp(B s — > K°K°) is very similar in both QCDF and pQCD (to NLO) approaches. It has 
been argued that the decay B s — > K°K° is a very promising place to look for effects of New Physics through 



the measurement of its direct CP violation fl57l.l58ll . For example, it was shown in (57j] that A CP {B S -> K°K°) 



which is not more than 1% in the SM, can be 10 times larger in the presence of SUSY while its rate remains 
unaffected. 



3 An early theoretical estimate yielded Z${B S —> K + K ) = (35 ± 7) x 10 6 using the measured B° — > K + K branching 
fraction l56ll . Based on QCDF and a combination of {/-spin and isospin arguments, a result of (20 ± 8 ±4± 2) x 



10 6 was obtained in 171]. 



17 



3. B^tjOtjC) 

The penguin-dominated rj^rj^ modes have sizable rates, especially B s —> rj'rj 1 ', the analog of B — > Krj' 
in the B s sector, has the largest branching fraction of order ~ 50 x 10~ 6 in two-body hadronic decays of the 
B s meson. The QCDF predictions in |0] within the S4 scenario are much bigger, 78 x 10~ 6 and 66 x 10~ 6 
respectively for r\r}' and r\'r]' modes. This is because Eq. ( 13.91 ) rather than (13.11b is employed there for 
describing the B s — > TjO transition form factors. One of us (CKC) found that the B s — > rj'rj' branching 



fraction can even reach the level of 1.0 x 10~ 4 in the residual final-state scattering model 13311 . It is evident 
from Table UTTI that the pQCD approach to lowest order predicts much smaller r\('7\0 rates even though 
the form factor Fq s71s (0) = 0.30 is used there. A recent pQCD calculation involving some NLO corrections 



from vertex corrections, quark loops and chormo-magnetic penguins exhibits some improvements [12]: the 
branching fractions of 7] 7], 7777' and rj'rj are enhanced from 8.0, 21.0 and 14.0 (in units of 10~ 6 ) to 10.0, 34.9 
and 25.2, respectively. The gap between pQCD and QCDF is thus improved. However, the NLO corrections 
calculated so far in pQCD are still not the complete results as some other pieces of NLO corrections such as 
hard spectator and annihilation have not been considered. It is important for the pQCD community to carry 
out the complete NLO calculations. 

Since the decays B s — ► rj 77 ( ) are penguin dominated and their tree amplitudes are color suppressed, 
their direct CP asymmetries are not large. 

4. B s — > KK 

The decays B s — ► %% proceed only through annihilation with the amplitudes 

A B s ^n+n- ~ V^Aj^bO^ « 2B n% b% (4.16) 
The predicted 3§{B S — > n + n~) = 2.6 x 10~ 7 in QCDF is consistent with the current upper limit of 1.2 x 10~ 6 



|29|]. Note that in the absence of power collections i.e. Pa = 0, the branching ratio will become too small, 



of order 5x10 



5. B s — > 71° Tj^ 

Since the isospin of the final state is / = 1 , the electroweak penguin is the only loop contribution that can 
contribute to the decays B s — > 7i°rj('\ in analog to the B — ► 7l~7t° transition. However, unlike the latter, the 
electroweak penguin amplitude in the former gains a CKM enhancement Xc /^i • Indeed, Pew dominates 
over C in B s — > 7l°Tj^ decays. It is well known that CP asymmetry of B — > k~k° is very small, of order 
10~ 3 . This is ascribed to the fact that the electroweak penguin there is very suppressed with respect to the 
color-suppressed tree amplitude C. On the contrary, CP violation of B s — ► 7F°7]C) is very sizable due to the 
dominant Pew- From Tables UTTI and II VI we see that the approaches of QCDF and pQCD have similar results 
for the rates of B s — > 7t°r]^ but quite different predictions for Acp(B s —>■ 7t°T]'). 
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TABLE V: Direct CP asymmetries (in %) in B s 
branching fractions and CP asymmetries for B d 



> PP decays via U -spin symmetry. Theoretical results of 
PP are taken from II17n . 



Modes 



#(10" 



Acp(%) 



Modes 



A CP {%)(U-spm) A CP (%)(QCDF) 



U d ~ 




IQ a+7.9+8.2 
Ly - J -4.&-6.2 


7 4+1.7+4.3 
'•^-1.5-4.8 




+ K + 7Z~ 


25.9 


90 7+5.0+3.9 
zu -' -3.0-8.8 


B d ~ 


+ n + n~ 


7 n +0.4+0.7 
' - u -0.7-0.7 


, 7 n+ 1.3+4.3 
1 '- u -1.2-8.7 


t~ 




-4.5 


7 7+1.6+4.0 
'•'-1.2-5.1 


t~ 




o £+3.8+3.8 
°'°-2.2-2.9 


-10 6+ 2 7 + 5 - 6 

iU '°-3.8-4.3 


t~ 


■,K°n° 


51.5 


T,/; -1 + 17.4+26.6 
JO - J -18.2-24.3 


t~ 




9 1+1.0+0.8 
Z - L -0.6-0.6 


_ 10 0+0.7+1.0 
iU - U -0.7-1.9 


t- 




0.77 


O 4O+0 4 +0 10 
u - w -0.04-0.04 


t- 


-+K+K- 


n 1n +0.03+0.03 
u - lu -0.02-0.03 





t~ 


+ % + %~ 









6. Test of U -spin and SU(3) symmetries 

There are five U -spin relations shown in Eqs. (14.51 ) and (14.71) . We have pointed out before that the relation 
(14.51 ) is experimentally verified. For other relations, we are still lack of the measurements of CP asymme- 
tries. Nevertheless, since the £/-spin and SU(3) symmetry breaking is already included in QCDF calcula- 
tions, we can test quantitatively how good the symmetry is. In Table [V] we show some of direct CP asym- 
metries in B s decays evaluated using the U -spin relations Eqs. ( 14.51 ) and ( 14.71 ) and theoretical inputs for the 
branching fractions of B d jS — ► PP decays and CP asymmetries of B d — > PP. We see that in general Aqp 
obtained by U -spin symmetry is consistent with that obtained from direct QCDF calculations. In jllll two 
parameters 

\A(B s ^7i+K-)\ 2 -\A(B s ^7i-K+)\ 2 



A = 



\A(B S - 
A C p(B d 



x+K-)\ 2 + \A{B s - 

m{B s 



K+K- 



+ 



7l-K+)\ 2 ' 
+ K-K+) T(B d ) 



(4 17) 

A CP {B S ^ jc-K+) ' &(B d ^> K+K-) t(S s ) ' 
are defined to quantify the U -spin violation through the deviation of Rt, from — 1 and A from 0. However, it 
is not suitable for the U -spin pair (B s — > K°K°, B d — > K Q K°) for which we find A « — 12. In this case, it is 
better to compare Acp(B s — > K°K°) obtained from the U -spin relation with the QCDF prediction as we have 
done in Table IVl 

As for the test of SU(3) symmetry, the first three relations in ( 14.121 ) are experimentally satisfied: 

5.0± 1.1 =5.16 ±0.22, 0.39±0.17 = 0.38±0.06, 24.4 ± 4.8 = 19.4 ±0.6, (4.18) 

where the branching fractions are in units of 10~ 6 and the data are taken from 14811 . For the last three 
relations of (14.121) we have 

-0.098±g;gjf , 26.1+}^ = 19.4±0.6, 0.004lg;^ = 0.009 ±0.025, (4.19) 



077+ 043 

u - u// -0.052 



where we have used the theoretical inputs for B s decays and experimental inputs for B d ones. Again, it 
appears that SU(3) symmetry relations are satisfactorily respected. 

7. Mixing -induced CP asymmetry 



Measurements of time-dependent CP asymmetries in neutral B s meson decays into a final CP eigenstate 
/ that is common to B s and B s will provide the information on two interesting quantities: mixing-induced 
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TABLE VI: Same as Table Hill except for the mixing-induced CP asymmetries Sf in B s —> PP decays. The 
parameter r\f = 1 except for K s (n° ,7] ,7]') modes where r/j = —1. 



Modes 


Class 


QCDF (this work) 


pQCD (LO) 


pQCD (NLO) 


SCET 


t~ 
t~ 


+ K S T1 


C 

c 


n n o+0.29+0.23 
U - UO -0.27-0.26 
rv 9fl +0.33+0.21 
U - ZO -0.44-0.30 


-0 61+ ' 24 

U - Di -0.20 

-0.43+°;§| 


-0.41+0-09 
-0-18^ 


-0.16 ±0.41 ±0.33 ±0.17 
0.82±0.32±0.11±0.04 
0.63±0.61±0.16±0.08 


t- 


+ K s ri' 


c 


n ns +0.21+0.20 
U - UO -0.17-0.16 


-0.68+°$ 


-0 46+ 012 

u - TO -0.23 


0.38±0.08±0.10±0.04 
0.24 ±0.09 ±0.15 ±0.05 


t- 
t- 
t- 


-> r/r/ 


p 
p 
p 


n 99+0.04+0.05 
u - zz -0.05-0.03 

004+ + 002 

n n7 +o.03+o.04 

u ' u/ -0.06-0.05 


0.28+°$ 

0.04 

03+ 001 
u - u:, -o.oi 


22+ a ° 4 

u - zz -0.03 

04+ 00 
u - u ^-o.oo 

02+ 000 
u - uz -o.oo 


0.19±0.04±0.04±0.01 

-0.026 ± 0.040 ± 0.030 ± 0.014 
-0.077 ± 0.061 ± 0.022 ± 0.026 


t- 


-> r/r/' 


p 


n m +0.00-0.00 
U - Ui -0.01-0.00 


04+ 00 


04+ 000 


0.041 ± 0.004 ± 0.002 ± 0.051 
0.015 ± 0.010 ± 0.008 ± 0.069 


t- 


->r\'r\' 


p 


rv ru+0.01+0.01 
O-Ot_0.01-0.01 


04+ 01 
u - u ^-o.oi 


05+ 000 
u - u;, -o.oi 


0.049 ± 0.005 ± 0.005 ± 0.031 
0.05 1 ± 0.009 ± 0.017 ± 0.039 


t- 


-> 7T TJ 


Pew 


O- /0 -0.23-0.47 


rv 1 t+0 ll 


rv <-vo+0 05 


0.45 ± 0. 14 ± 0.42 ± 0.30 
0.38 ±0.20 ±0.42 ±0.37 


t- 




Pew 


n oe+0.03+0.04 
O-°°-0.15-0.29 


-0 17+ 008 

u - i '-0.09 


18+ 012 


0.45 ±0.14 ±0.42 ±0.30 
0.38±0.20±0.42±0.37 


t- 
t- 




ann 
ann 


n 1 <r +0.00+0 
u - iJ -0.00-0 

n 1 <r+o.oo+o 

u - iJ -0.00-0 


14+ 012 
14+°- 12 

U ' i ^-0.06 


09+ 02 

rv ns+ 00 
o-o«_o.OO 





CP asymmetry Sf and direct CP violation A f which can be expressed as 

1-IA/l 2 21mA/ 
1 + N 2 ' 1±|A/| 2 ' 

q Bs A(B s ^f) _ V t * b V ts A(B s ^f) 



where 



(4.20) 



(4.21) 



f p Bs A(B s ^f) V tb V t *A(B s ^fY 

Now let q B ^ / p B = e 2 '^ s and 

A(B s ^f) = A l e i ^ Al+s ^+A 2 e i ^ A2+S2 \ 
A(B S /) = T] f (A l e i{ -' t ' Ai+Sl) +A 2 e i ^ A2+&l) 

where CP\f) = f]f\f) with r// = 1 (— 1) for final CP-even (odd) states, 0aia2 are weak phases and 5^2 
strong phases. It follows that (see e.g. 159(1 ) 

1 + re K<l>i-<h) e iS 



(4.22) 



1 u i + re -i(<h-<h) e lS 
with 0i 2 = 0ai A2 + ft, S = 82 - -1 and r = A 2 /A l . 



(4.23) 
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For B s decays, the phase f5 s due to the B s — B s mixing is very small in the SM, of order 1°. For the 
decays B s — ► K°K°, T]T], 77 ry', T^'t]' dominated by penguin diagrams (tree contributions to TJ^TjO are color 
suppressed), r ~ and the phase 0ai due to V^V^* or V t bV t * is also very small. Consequently, 5/ are small for 
penguin-dominated B s — > PP decays except for B s — > K + K which receives a tree contribution with 0^2 = 
We see from Table|Vl]that QCDF, pQCDF and SCET all predict S S ^ K+K - 0.20. Recently, both CDF 
and DO J2I] have reported fits to angular and time distributions of flavor-tagged B s — ► J/yty decays which 
favor a larger value of f} s deviated from the SM by l-2a effects. If this is the case, then mixing-induced 
CP violation in B s — ► K K G ,r\r] ,r\r}' ,r\' r\' could be sizable. Hence, these modes offer rich possibilities of 
testing New Physics beyond the SM. 

Due to the large magnitude and strong phase of a-i induced from soft power corrections to the color- 
suppressed tree amplitude, for example, ax(K%) = 0.77e~' 52 ° (or 0.41e~' n ° before corrections), 4 we find 
that such corrections will flip the sign of Sf into the positive one for the color-suppressed decays B s — > 
7^5(71°, 77,77'), while they are all negative in the pQCD approach. Recently, it has been claimed that in the 
pQCD approach there exist uncanceled soft divergences in the kj factorization for the nonfactorizable B 
meson decay amplitudes 16011 . This will enhance the nonfactorizable color-suppressed tree amplitudes. It 
remains to check if the signs of Sg s ^ Ks (-p ^ in pQCD will be flipped again under this "#2" enhancement. 



V. B s -> VP DECAYS 
A. Branching fractions 

The tree-dominated decays B s — > K* + 7i~ and p K + with the amplitudes 

A(B S ~^ K* + 7t~) « A K *„(S pu a- l +a 4 -r*a 6 ), 

A(8 s ->p-K + ) Rs A Kp (8 pu a l +a 4 + r P x a 6 ), (5.1) 

have branching fractions of order 10~ 5 . Since A K , n = » f n A% sK * (OJmJ, and A Kp = X^ K P"> » 

fp F sK (®) m B t see Ec l- <l231) 1. it is clear that the p K + mode has a rate larger than K* + n~ due to the 
hierarchy of the decay constants f p S> f % . The penguin-dominated B s — > VP decays such as K* K + and 
^*o^o ^ ave rates sma n er than the counterparts in the PP sector as the amplitudes are proportional to a 4 + 
r^T £?6 or CI4 — r^as for the former and a 4 + r^a(, for the latter. Since a\ and a(, are of the same sign and 

> , it is evident that the interference of the penguin terms is constructive for PP and either destructive 
or less constructive for VP. The decay B s — > (j)K° is dominated by the b — > d penguin transition and its rate 
is thus much smaller compared to b — > s dominated B s — ► K*K decays. 

We see from Table I VIll that the pQCD predictions for the color-suppressed tree-dominated decays B s — > 
K*°7l ,p°K , (oK°,K*°ri' are one order of magnitude smaller than QCDF and SCET in rates. For example, 
3B(B S — > p°K°) is predicted to be of order 1.9 x 10~ 6 in the approach of QCDF, but it is only about 0.08 x 
10~ 6 in pQCD. The calculated branching fractions in pQCD for K*°ri and some of the penguin-dominated 
decays e.g. B s — > K* + K ,K*°K° ,<pK° ,§r\' are also much smaller than QCDF. In the following we will 



4 In the B u c i systems, Q2{Kk) = 0.5 le ,58 ° (or 0.27e ,17 ° before corrections). 
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comment on the decays B s — ► 07] 0. While the QCDF approach leads to ^(S. s — ► 07]') > ^(B v — > 0?]), 
pQCD and SCET predict very different patterns: ^(B, — ► 07] ) > ^(fi s — >• 07]') in the pQCD approach and 
it is the other way around in SCET (see Table IV"iTl ). We should stress that the decay rate of B s — ► 07]' is 
sensitive to the form factor Aq s ^(0). The decay amplitudes of B s — > 07] ^ are given by 

A(fl s ->07]) = cos0A(5 4 . -> 07] 9 ) -sineA(5, ; -> 07],), 

A(S,^07] / ) = sin0A(B,^07] ? )+cos0A(B i ^07],), (5.2) 

with 

A(S,^07]. S ) = A ^< + <)+^K + <), 

\/2A(S s -> 07] ? ) = A^ 9 (S pu a 2 + 2a 3 p ). (5.3) 

Since a^((j)rj s ) = a$ — r^a^ and a|(7] 5 0) = «4 + are of opposite sign (numerically, a|(07] s ) « 0.038 
and a|(7] s 0) ~ —0.033), there is a cancelation between the two penguin amplitudes of fi s — > 07] v . Note that 
a|(07] v ) and (tj s 0) also are of opposite sign. It turns out that the sign of A(B S — > 07],) depends on the 
form factor A () (0). For A (0) = 0.32 as employed in the present work, B s — ► 07]. v and B s — > 07]^ will 
contribute constructively to — > 07}' so that — > 07]') = 2.2 x 10~ 6 . However, if we use the sum-rule 
prediction (0) = 0.474 from Eq. (l33i then a near cancelation between B s — > 07] s and B. v — > 07]^ occurs 
in the decays S s — ► 07]', so that its branching fraction, of order 10 -7 , becomes very small. Hence, it is very 
important to measure the branching fractions of B s — ► 07] ^ to gain the information on the form factor A^ . 

One unique feature of the B s decays is that there exist several modes dominated by electroweak penguins: 
B s — ► 7r°7]^, 071°, p ?]'') and 0p°. The isospin for the final states of these decays is / = 1 and hence the 
electroweak penguin is the only loop contribution that one can have. It dominates over the color-suppressed 
tree contribution due to the large CKM matrix element associated with the electroweak penguin amplitude. 
Since a large complex electroweak penguin amplitude due to New Physics is also a possible solution to the 
B — > Kk CP puzzle, it has been advocated that this hypothesis can be tested in the decays B s — ► 07T°, 0p° 



whose rates may get an enhancement by an order of magnitude 16111 . 



B. Direct CP asymmetries 

Direct CP asymmetries of B s — > VP decays estimated in various approaches are summarized in Table 
IVIIII In QCDF calculations, the signs of CP asymmetries for color-suppressed tree-dominated decays B s —> 
K*°7l ,p°K , (oK° and K*°ri are governed by the soft corrections to a 2 [see Eq. (13.26H . We see that QCDF 
and pQCD results agree with each other in signs, whereas SCET predicts opposite signs for these modes. 
Since the corresponding rates of these decays are very small in pQCD, as a consequence, the CP-violating 
asymmetries predicted by pQCD are very large, of order 0.50 or even bigger. 

In the pQCD approach, the penguin-dominated decays B s — > K°(j),K* K ,K* K have no direct 
CP asymmetry as their decay amplitudes are governed by one type of CKM matrix elements, e.g. V t bV* d for 
the first mode and V t bV t * for the last two. As noticed before for the decay B s — ► K°K°, NLO corrections from 
penguin loop interactions can bring a weak phase necessary for a non-vanishing CP violation. Therefore, 
it is important to carry out pQCD calculations to NLO for those three modes. In the approach of SCET, 
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TABLE VII: CP-averaged branching fractions (in units of 10~ 6 ) of B s — ► PV decays calculated in various 
approaches. The LO pQCD predictions are taken from jllll . while two different sets of SCET results are 
quoted from 11411 . 



Modes 


Class 


QCDF (this work) 


PQCD 


SCET 1 


SCET 2 






r 


7 n+0.4+0.5 
•— 0.7— 0.7 


7 .-+2.9+0.4+0.5 
• — 2.2— 0.5— 0.3 


c q+0.5+0.5 
J - -0.5-0.5 


fi fi +0.2+0.7 
u-u — 0.1— 0.7 


B° s - 


^p-K+ 


r 


l4.1 + \i+°A 

-1.9—1.3 


17 S+7-7+1.3+1.1 
•°— 5. 6-1.6-0. 9 


7 fi +0.3+0.8 
' — 0.1— 0.8 


in 9+0.4+0.9 
—0.5—0.9 


B°s- 


->K*°n° 


c 


rv oq+0.80+0.84 
"■° -0.34-0.35 


rv r, 7 +0.02+0.04+0.01 
u ' u -0.01-0.02-0.01 


rv qrv+0.07+0.10 
"•'"— 0.01— 0. 11 


i rv 7 +0.16+0.10 
u -0.15-0.09 


B° s - 


^p°K Q 


c 


i q+2.9+1.4 
-"-0.9-0.6 


o ns+ 002 + 007 + 001 

U.UO_q 02-0 .03-0 .00 


2 rv+0.2+0.2 
-0.2—0.2 


rv o 1+0.05+0.08 
• -0.02-0.09 


B° s - 


^(0K° 


c 


i fi +2.2+1.0 
-0.7-0.5 


n i ^+0.05+0.07+0.02 
■ -0.04-0.03-0.01 


0.90+^+Oj? 

—0.01—0.1 1 


i n+0.1+0.1 
-0.1-0.1 


B° s - 


->K*-K+ 


p 


10 3+^+* ~ 


(r +1.7+1.7+0.7 
"•— 1.5— 1.2-0.3 


o 4 +4.4+1.6 
" — 3.4— 1.3 


q c+3.2+1.2 
?-J -2. 8-1.1 


B°s- 


-> K*+K~ 


p 


11.3+™+®! 

— 3.5—5.1 


4 7 +l.l+2.5+0.0 
-0.8-1.4-0.0 


—3.7—1.4 


— 3.2—1.2 


«?- 


-> ^*°^° 


p 


10.5+3M.1 

— 2.8—4.5 


7 1+2.5+2.1+0.0 

/ • J 1 *~7 1 O A A 

— 1.7—1.3—0.0 


7 q+4.4+1.6 
■ —3.4— 1.3 


q n+3.2+1.2 
—2.8—1.0 


B° s - 


-> tf*°£ 


p 


— 3.6—4.8 


A ^+0.7+2.2+0.0 
—0. /— 1 .4—0.0 


o 7 +4.4+1.6 
— 3.5—1.4 


q 4 +3.7+1.4 
— 3.1 — 1.2 


B° s - 


^(j)K° 


p 


0.6+^+n 

—0.2—0.3 


n 1^+0.04+0.09+0.02 

\J ■ -L V_J A AO A A ,4 A A 1 

— 0.03—0.04—0.01 


rv 44+O.23+O.O8 

V / • I I A 1 O A A^ 

— O.lo— 0.07 


n SzL +0.21+0.08 
—0.17—0.07 


8}- 




Pew 


n i O+0.02+0.04 
— 0.01 — 0.02 


n 1^+0.06+0.02+0.00 

\J ■ -L VJ A A C A A^ AAA 

—0.05—0.02—0.00 


n n7 +0.00+0.01 
—0.00—0.01 


n m +0.00+0.01 

V/.V/Z7 A A A A A1 

—0.00—0.01 


B° s - 


-> p + 7r~ 


ann 


n rp+O.OO+O.Ol 
—0.00—0.01 


rv 99 +0.05+0.04+0.00 
—0.05—0.06—0.01 








-> p~7T + 


ann 


n m +0.00+0.01 
—0.00—0.01 


n 9^+0.05+0.05+0.00 

\_/ . Z*^T A AC f\ f\/Z A A 1 

—0.05 — 0.06—0.01 






B°s- 




ann 


rv 2+ 00 + 01 
—0.00—0.01 


r. 9-3+0.05+0.05+0.00 
-0.05-0.06-0.01 






B°s~ 


->^*°r/ 


C 


^+0.33+0.35 
—0.14—0.17 


n i 7 +0.04+0.10+0.03 
— 0.04—0.06—0.01 


i 7+0.3+0.2 
■ —0.3—0.1 


n ^ 9 +0.14+0.07 

U.UZ,_a |4_a no 


B}- 




C 


^•^-0.30-0.41 


r,Q+ 02 + 03 + 01 
^•^-0.02-0.02-0.01 


u - w -0.26-0.11 


n 07+O.35+O.IO 
u - 0/ -0.32-0.08 


B°s- 


+ 071 


P 


i n +l. 3+3.0 
1 - u -0.1-1.2 


-3 -r + 1. 5+0.8+0.0 
1.0-0.6-0.0 


n <rq+2.02+0.12 
u - J ^-0.59-0.10 


Q4+1.89+0.16 
"•^-0.97-0.13 






P 


9 9+4.5+8.3 
z - z - 1.9-2.5 


n iQ+0.06+0. 19+0.00 
u - i ^-0.01-0.13-0.00 


7 n+7.7+1.6 
'• J -5.4-1.3 


4 O+5.2+0.7 
^•-'-3.6-0.6 




-> 0)7] 


P,C 


n r.i+0.12+0.06 
U - UJ -0.02-0.01 


n nA+° °3+o o5+o.oo 

^•^-0.01-0.02-0.00 


rv r,4+0.04+0.00 
u - UH -0.02-0.00 


rv r.r, 7 +0.01 1+0.001 
u - uu/ -0.002-0.001 




-> air/' 


P,C 


n ,^+0.27+0.15 
u - iJ -0.08-0.06 


n AA +0. 18+0. 15+0.00 
^•^-0.13-0.14-0.01 


rv rvrvi +0.095+0.000 
u ' uui -0.000-0.000 


r, 9 n +0.34+0.02 
u - zu -0.17-0.02 


B° s ~ 


^p°r> 


Pew 


n 1n +0.02+0.02 
u - iu -0.01-0.01 


n r.A+ 03 + 01 + 00 

^"-oxg-o.oi-o.oo 


n ns+ 04 + 01 
^•^-o.os-o.oi 


n nA+ 03 + 00 

L, - L,o -0.02-0.00 


B°s- 


+ pY 


Pew 


rv i/:+0.06+0.03 
u - iD -0.02-0.03 


rv i O+0.06+0.02+0.00 
u - i:> -0.04-0.02-0.01 


n nn^+ 082 + 000 

U.LIUJ_ 000 _ .000 


rv 14+0.24+0.01 
u - i ^'-0.11-0.01 



CP asymmetries of the decays B s — ► and B A . — > p°(ft))(T] , T]') also vanish. As explained in Il4ll . there is 
no charming penguins in these 5 channels and hence no direct CP violation due to the lack of strong phases. 

We use this chance to clarify one misconception about CP violation under isospin symmetry. The isospin 
of the final-state is / = 1 for B s ^071°, p°rj^ and 1 = for B s — ► (0, ft))T]( ; . One may argue that there is no 
CP violation for these decays as they have only one isospin strong phase (see e.g. 06211 ). 5 On the contrary, 



we found large direct CP-violating effects in some of above decays (see Table IVIII1) . The point is that 
isospin phases should not be confused with other possible strong phases in each of topological amplitudes. 
In our study, CP asymmetries of B s — > pTjO are large since the electroweak penguins dominate over the 
color-suppressed tree amplitudes. 



5 By the same token, it has been (wrongly) claimed that the direct CP asymmetry is strictly zero in the charged 
B — ► 7T _ 7T° decay. 
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TABLE VIII: Same as Table I VIll except for the direct CP asymmetries (in %) in B s — > decays. 



Modes 


Class 


QCDF (this work) 


PQCD 


SCET 1 


SCET 2 


tf- 




T 


-24.0+^+" 

— 1.5— 3.9 


10 n+2.5+2.7+0.9 

—2.6—3.4—1.4 


q q+17.2+0.9 
— 16.7—0.7 


_ 12 4+17.5+1.1 

J- Z. .— ic It 

— 15.3—1.2 


t- 




T 


— 2.1 — 1 1 .6 


i/l 9+2.4+2.3 + 1.2 
— 2.2—1.6—0.7 


1 1 0+17.5+1.2 

ll.O ")AA T T 

20.0 1. 1 


10 S+9.4+0.9 

lu.O 1AO 1 A 

— 10.2— 1 .0 


t- 


^K*°n° 


C 


ryr ^, + 10.8+42.2 
— 10.9— 36.7 


47 ,+7.4+35.5+2.9 
— 5. / — 29.5— /.0 


99 Q+33. 1+2.1 

^■Z>.7 a t\ i in 

—40.2—1.9 


10 A +18.6+0.8 
1 -J . 1 100 it 

— lo.o— 1 .2 


t- 


^p°K Q 


C 


90 q+14.6+25.0 

Z»0.7 1 A C TO *1 

— 14.5—23.7 


7^ /1+6.4+16.2+2.2 

1 J 'i 1 1 ^ ai i a 

— 11. /— 4/.5— 3.9 


19 n +2o.i+i.o 

1 ^ ■ W 1 f\ f f\ ^7 

— 19.6—0.7 


09 ^+30.7+2.7 
—23.4—2.9 


t- 


-> coK Q 


C 


09 n +18.9+23.6 
— 17.5 — 26.2 


c 9 1+3.2+22.7+3.2 

*J Z- .1 f\ A 1CT *> / \ 

—0.0—15.1—2.0 


?4 A+33.7+2.2 
—41.4—2.0 


10 9 + 16.4+1.2 
— 17.0—1.7 


t- 


^K*-R+ 


P 


1 1 n +0.5+14.0 

1 1 . VJ f\ A IOO 

—0.4— 15.5 


-36.6+"+2-8+|-3 

—2.3—3.5—1.2 


1 1 9 + 19.1 + 1.3 
— 16.2— 1 .3 


19 o+H.4+0.8 

l^--5 11T AO 

— 1 1.3— 0.O 


t- 




P 


9S <r+9.2+16.3 

Z.-J .>J OO 111 
— 5.5— 1 1 .3 


cc ^+4.4+8.5+5.1 
—4.9—9.8— 2.5 


7 1+11.2+0.7 
— 12.4—0. / 


g fi +13.0+0.7 
— 13.5—0.9 


t- 


+ T°K° 


P 


n /1Q+0.08+0.09 
— U.U / — U. Iz 











t- 


-> /r°}f 


P 


n 1n +0.08+0.05 

\J • 1-\J t\ c\"7 a m 
—0.0/ — 0.02 











t- 




P 


n 9 + I.2+O.6 

J .Z, 1 A IO 

— 1 .4— 1 .j 





n n +5.3+0.3 

—4. / — U. j 


9 9+3.O+O.I 
— z.y— u. 1 


t- 


■*(j>n° 


Pew 


09 9+10.9+ 9.0 

— 14.0—55.3 


1 , ,+0.3+2.1+1.5 

I J. J A/1 1 "7 A. "7 

—0.4— 1. 1 — 0. / 








t- 


-> p + 7T" 


ann 


in 9+0.8+12.7 

— 0. /— lZ.5 


A -r +0.0+2.9+0.6 
—0.6—3.5—0.3 






t- 


-> p~7T + 


ann 


. 1 1+0.7+13.9 

11.1 AO ic i 

— 0.5— 15. / 


1 -3+0.9+2.8+0.1 
—0.4—3.5—0.2 






t- 


4 P V 


ann 





1 7+0.2+2.8+0.2 
1 * ' no 1 ^ n 1 

— U.5 — J.D — U . 1 






t- 




C 


— 19.2—64.5 


ci 7+6.2+14.1+2.0 
— 6.4—12.4—3.3 


9<r 7+23.4+2.0 
—22.0—1.3 


-T9 7+28.1+2.6 

— 22.5—3.9 


t- 




C 


fry ,-+6.0+24.7 

— 5.5—20.2 


ci 1+4.6+15.0+3.2 
—6.6— l5.z— 4. 1 


9+63.3+3.1 
—49.4— J. 


09 1 +22.8+2.6 

J^.l no i 1-7 

— 23.2— 1 . / 


t- 




P 


19 A +14.1+64.9 

— 5.7—39.8 


1 0+O.O+O.6+O.I 

l.O /AT / \ / r\ 

—0.1 — 0.6—0.2 


91 n+53.5+2.5 
— 83.2—2.6 


ifi q+13.8+1.6 

1 VJ .7 1 t i 

— 18.3 — 1.6 


t- 


->017' 


P 


it q+15.4+28.5 
4.2-89.7 


7 O + 1.5+1.2+0.1 
'•°-0.5-8.6-0.4 


A 4+5. 3+0.6 
^•^-7.1-0.6 


7 O+5.0+0.8 
' -0-4.9-0.8 


t- 


-> COT] 


PC 


fiA Q+ 24 - 4 + 44 - 


\f. 7+5.8+15.4+0.8 
iu - '-3.2-19.1-1.7 








t- 


-> COT]' 


PC 


/t+4.4+10.4 
Ji '- H '-3.0-11.7 


7 7+0.4+4.5+9.4 
'•'-0.1-4.2-0.4 








t- 




Pew 


7 <r 7+15.3 + 13.3 
/J - '-17.6-37.5 


q 7 + 1.0+2.8+0.4 
^■ z -0.4-2.7-0.7 








t- 


^p°r/' 


Pew 


07 4 + 3.4+ 5.7 
°'- H -10.6-30.3 


9<r + 1.3+2.8+3.4 
z,j.o_ 2 _ 3 6 _j 5 









C. Test of U -spin and SU(3) symmetries 



The pairs related by £/-spin symmetry are (63|]: (B d — > K*~n + , B s ->• p~K + ), (B d ->• ^ _ p + , -> 
(5d P"^ + 5 5 -v ^*^ + ), (S d P+7T-, 5, /?*+/?-), (S £/ B s -» p-7T+), (S d 

/?*+/?-, B, -> p+7T"), (Bj -f /?*°/?°, B 5 -> *T* Z°), -> /?*°/?°, -> /?*°/?°). Note that unlike and 

modes, B 5 -> K*°K° and B s -> K°p° are not related to B d -> p°K° and -> K*°K°, respectively. Direct 
CP asymmetries of the pairs listed above are related by U -spin symmetry in analogue to Eq. (14.51 ) or Eq. 
(I4.7I ). The test of £/-spin symmetry in B s — > VP decays is shown in Table llXl It turns out that £/-spin 
symmetry is in general acceptable. 

Just as B s — ► PP decays, under the approximation of negligible annihilation contributions to tree- 
dominated decays and kee ping on ly the dominant penguin annihilation terms in penguin-dominated decays, 



SU(3) symmetry leads to ||46U49ll 



A(B S K* + n-) ^A(B d P+7T-), A(B S p-/?+) » A{B d -> p-jr+), 

A (5, -> /T+*r) «A(B rf ^p+^), A(S^/T*-/r + ) »A(S d ->Z*-7T + ). (5.4) 
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TABLE IX: Direct CP asymmetries (in %) in B s —> VP decays via £/-spin symmetry. 



Modes 



^(i(r 6 ) 



A CP {%) 



Modes 



A CP (%)(U-spm) A CP (%)(QCDF) 



B d 
B d 
B d 

B d 

b2 



JT+lsr 

K*-R+ 

K*°K° 

K*°K° 

p + n~ 

p~n + 



q 9 +1.0+3.7 
^■ z -1.0-3.3 



+5.7+7.4 
-2.8-4.5 
+0.01+0.02 
-0.01-0.02 
+0.01+0.04 
-0.01-0.03 



n 7f) +0.18+0.28 
u -' u -0.15-0.25 



8.6 

0.08 
0.07 



0.47 



hO.36+0.43 
-0.17-0.27 



q 9 +0.4+0.5 
-0.7-0.7 
+ 1.1+0.9 
-1.5-1.1 



-12.1 
31.9 



0.5+12.6 
0.5-16.0 
+ 11.5+19.6 
-11.0-12.7 



a 7 +0.1+4.7 
^•'-0.2-2.7 
c s +0.2+7.0 
J - J -0.2-5.5 
it c + 1. 6+1.4 
iJ - J -1.7-2.3 



15.9 



-3.5 
-22.7 
4.4 



-1.3- 
-1.7- 
-0.9- 
-1.1- 
-0.3- 



-0.3-6.8 



t 
t 



p~K+ 
■K*+n- 

■ p + 7T _ 

■ p~K + 

■ K*°K° 
■K*°K° 

K*+R- 
K*K+ 



9.1 

-39.6 
18.9 

-20.5 
0.86 
0.17 
19.0 

-6.6 



1 1 7+ 35 
n./_ 2 .i 



10.1 
11.6 

-24 +L2+7 - 7 

ZH - u -1.5-3.9 

in 9+0.8+12.7 
lu ' z -0.7-12.8 
_i i 1+0.7+13.9 
ii - i -0.8-15.7 
n .iq+0.08+0.09 
u '^-0.07-0.12 
n 1n +0.08+0.05 
u - iu -0.07-0.02 
9c r <r+9.2+16.3 
ZJ ' J -8.8-11.3 
-0.5 + 14.0 
-0.4-18.8 



-11.0 



Thus, we have the relations 

^- — w / __j,_L 



and 



3&{B S 

A C p{B s - 
A C p{B s - 



>K* + 7t-); 
> K* + K 



iA CP (B d - 
*Acp(B d - 



>p + K ), 



>p + K ), 

> P +n, 



m(B s - 

A CP (B S - 
A C p(B s 



* K*~K + ) t 



• @{Bd ■ 



p 7T+), 



iA CP {B d - 
~A CP {B d 



■p 7T+), 



Numerically, 



14.7^ 
10.3+^ 



-q 0+3.8 



7 c+0.6 j_ q 9+0.6 
'•"-1.0 — ^- z -1.0 ' 

11.3+^ = 8.6+^, 
for branching fractions in units of 10~ 6 and 

- 24.0+1* = -22.7+1 j, 

9C c+18.7 j_ o, q+22.7 

ZJ ' J - 14.3 ~~ Ji-y_ 16.8 ' -lo.o -io.u' 

for direct CP asymmetries in %. Hence, the above SU(3) relations are generally respected. 



i i 7 +10.7 j_ 4 4 +5.8 

_ii n+14.0 j_ _ ]2 1+12.6 
11 - u -18.8 ~~ 1Z ' 1 -16.0 ' 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



D. Mixing-induced CP asymmetry 

As discussed before, due to the tiny phase in the B s — B s mixing and in the CKM matrix element V c bV* s 
or VtbV* s , mixing-induced CP violation Sf is expected to be very small in the penguin-dominated B s . — > 0t?' 
decays. This is indeed borne out in all model calculations. The b —> dg penguin-dominated decay B s — > Ks<p 
has a large mixing-induced CP asymmetry due to the fact that the CKM matrix element V u bV* d has a weak 
phase — y. More specifically, 

A(B S -> K s $) - VubKi I A ^«3 + A^ K {< + ft")] + V cb V* d [A^ a L 3 +A^ K {a c 4 + j3 3 c )] . (5.9) 
To the approximation that a^ A « a 3 % and j6 3 c « j3 3 ", it is clear that MB S -> ^ 5 0) oc v, fe V^ = |V r ffc V r ^|e' /3 . Thus, 



5 b s ^4> ~ - sin2 (A +P) = -°- 71 for ft ~ !° and P = 21.58° |37y. In the pQCD approach, this decay is 
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TABLE X: Same as Table IVDl except for mixing-induced CP asymmetries Sf in B s — > PV decays. The 
parameter r\f = 1 except for K s (p°, ft), 6) modes where r\f = — 1. Note that the error estimate of Sg s _> K a is 
not available in the pQCD calculation II 111 . 



Modes 


Class 


QCDF (this work) 




pQCD 




SCET 1 


SCET 2 


t~ 




C 


n 9Q+0.23+0.16 

-U.Z4-U.zl 


-0.571°$: 


-0.51+0.02 
-U.iy— u.Uj 


0.99; 


fO.00+0.00 
-0.03— O.Ol 


a ai+0.22+0.17 
—0.1 /— O.lz 


t~ 


^K s a> 


c 


a Q9+0.03+0.08 
— U.U7— 0.15 


-0.63+"-^- 
— u.uy- 


-0.28+0.01 
-O.l 1 — O.Oz 


-o.ii; 


f-0.28+0.18 
-O. zz— 0.14 


n qc+o.02+0.00 

— 0.04 — 0.01 


t~ 




p 


fiQ+0.01+0.01 

— o.oi — u. oi 




-0.72 


0.09; 


f0.04+0.01 

/i /n a A i 
-0.03—0.01 


a i -,+0.02+0.01 
—0.02— 0.01 


t~ 




p 


a 91+O.O8+O.6I 
u ' z -0.11-0.25 


-0 03 +0 02 ~ 
u - UJ -o.oi- 


-0.07+0.01 
-0.20-0.02 


-0.39; 


(-0.43+0.04 
-0.15-0.04 


a 9^+0.35+0.02 
u - ZJ -0.16-0.02 


t- 


+w 


p 


a aq+0.05+0.48 
U - UO -0.06-0.81 


a oo+ °- 


-0.02+0.00 
-0.02-0.00 


-0.07; 


1-0.06+0.01 
-0.06-0.01 


a ia+0.07+0.01 
u ' iu -0.05-0.01 


*- 


-> 0)7] 


P,C 


a 7fi +0.16+0.52 
u - /o -0.03-0.22 


u - uz -0.03- 


-0.02+0.00 
-0.08-0.00 


-0.62; 


1-0.41+0.08 
-0.18-0.12 


a Qo+0.04+0.03 
u -^ J -0.98-0.04 


t- 


-> 0)7]' 


P,C 


a 8/L +0.06+0.04 
u -°^-0.05-0.03 


u - 11 -0.00- 


-0.04+0.02 
-0.04-0.03 


-0.25; 


hl.23+0.10 
-0.74-0.16 


i aa+0.04+0.01 

1 - uu -o.oo-o.oo 


t- 




Pew 


a a a+0.04+0.32 
u - HU -0.10-0.53 


-0 07+ 01 - 
u - u/ -o.oi- 


-0.08+0.02 
-0.09-0.03 


0.89; 


1-0.00+0.04 
-0.00-0.05 


A QA + 0.00 + 0.02 

"•^-0.00-0.03 


t~ 




Pew 


a o<r+0.09+0.22 
u - ,J -0.16-0.40 


1 5 + 006 - 

u ' iJ -0.06- 


-0.14+0.01 
-0.16-0.01 


l.oo; 


1-0.00+0.00 
-0.06-0.01 


60 +a30+0 03 

u - uu -0.53-0.03 


t- 


-+p r/' 


Pew 


a AJ r+0.05+0.30 
u -^ J -0.13-0.35 


u - i0 -o.oo- 


-0.10+0.04 
-0.12-0.05 


0.95; 


1-0.00+0.02 
-1.60-0.02 


a a, +0.75+0.10 
V -^ L -0.75-0.15 


t~ 




ann 


a 6< r+0.03+0.00 
u - D - 5 -0.03-0.00 


1 q+0.00- 
u - iy -0.00- 


-0.02+0.01 
-0.02-0.02 









dominated by the (S — P) (S + P) penguin annihilation process with the CKM matrix element proportional to 
V t bV* d . Therefore, both QCDF and pQCD predict Sg^^ ~ ^(0.70). (However, no error estimate is done 
in the pQCD calculation 111 in .) On the contrary, the SCET result of Sj> ~ 0.09 or —0.13 is dramatically 
different from the QCDF and pQCD predictions. As explained in 11411 . charming penguin contributions 
to B s — > Ks<j) dominates over penguin operators and the CKM matrix element associated with charming 
penguins is V c bV* d . Hence, S^^ Kst p = — sm2j5 s = —0.03 is predicted by SCET when penguin contributions 
are neglected. It should be stressed that although both QCDF and pQCD approaches have similar results 
for Sjj s ^ Kst j,, they differ in the prediction of Acp(B s — > K$(j)): it is of order —0.03 in QCDF and vanishes in 
pQCD for reasons mentioned above. 

The study of CP violation for B s — > K* + K~ and K* K + is more complicated as K^K^ are not CP eigen- 
states. The time-dependent CP asymmetries are given by 

r(5j(f ) -* K*±KT) - V(B%t) -> K*±K^) 

r(flj(f) -» K*±K*) +T(B° s (t) -> K*±K^) 
= (5±A5)sin(Am s /)-(C±AC)cos(Am s /), (5.10) 

where Am s is the mass difference of the two neutral 5 S eigenstates, S is referred to as mixing-induced 
CP asymmetry and C is the direct CP asymmetry (C = —Acp), while AS and AC are CP-conserving quanti- 
ties. In writing the above equation we have neglected the effects of the width difference of the B s mesons. 
Defining 

A+_ = A(B° ^K* + K~) , A_ + =A(B° ^K*~K + ) , 

A- + = A$, ^ K*~ K + ) , A+_=A(flJ ^K* + K~), (5.11) 



and 



A + - = ^^, A_ + = ^^±, (5.12) 
P Bs A +- Pb s a -+ 
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TABLE XI: Various CP -violating parameters in the decays B° s — ► K^K^. SCET results are quoted from 

i3. 



Parameter 



QCDF (this work) 



c 
s 

AC 
AS 



n iQ+0.03+0.14 
u - ly -0.04-0.11 



0.04+0.15 
-0.04-0.14 
+0.01+0.13 
-0.01-0.09 



n O o +0.12+0.46 
-"■"J_0.14-0.49 



-0.08 

-0.05 



0.33 



-0.09+0.30 
-0.10-0.48 



SCET 1 



SCET 2 



0.10+0.00 
-0.11-0.00 

0.07+0.01 
-0.07-0.01 



a nq+0.1 1+0.01 
~ u - uy -0.10-0.01 



0.02 

-0.02 



0.01 
0.02 



0.38 



-0.07- 
-0.07- 



0.04 
0.04 



+0.09+0.00 
-0.09-0.00 
+0.05+0.01 
-0.05-0.00 



n i i +0.09+0.01 
u - 11 -0.09-0.01 
n A 1+0.05+0.03 
U - H1 -0.05-0.03 



we have 



and 



I 2 IA+ l 2 -|A+ I 2 1-IA-+I 2 |A_+| 2 -A +| 2 

C + AC=- — ' + ' ' + ' C-AC=- — l pr^\j= ' , ? = + ' (5.13) 
1 + |A+_| 2 |A + _| 2 + |A+_| 2 ' 1 + |A-+| 2 |A_ + | 2 + |A_+| 2 ' 

5 + AS" 2lmX+ - _ 2lm ( e2iPsA +- A *+-) 



1 + |A + _| 2 |A + _| 2 + |A + _| 2 ' 
2ImA_+ 2Im(e 2 *A_ + A* , ) 

s - ASs TW^? = |A-Vn-+P <5 ' 14) 

Hence we see that AS describes the strong phase difference between the amplitudes contributing to 
B° s -> K* ± K T and AC measures the asymmetry between F(B° S -» K* + R-) +T(B° S -> K*-R + ) and F(B° -> 
K*-R+) +T(lf s -> K*+K-). 

Next consider the time- and flavor-integrated charge asymmetry 

_, _ |A+-| 2 + |A + _| 2 -|A_ + | 2 -|A_ + | 2 

^ K = " '2,11. |2,|_ .12 ,11 .12' < 5 - 15 > 



|A + _| 2 + |A + _| 2 + |A_ + | 2 + |A_ + 



Then, following 13711 one can transform the experimentally motivated CP parameters s^k*k and Ck*k into 
the physically motivated choices 



k--+| 2 -1 |jc+-| 2 -1 



Ak*-k+ = (5-16) 



with 



Hence, 



1 k*+k- 



\k-+\ 2 + V tt n \k+-\ 2 + \ 

<hs_ ^± -+ _ A+- 

P* A +- ' Pb s A -+ 



r(_ff -> _;*+_: ) - r(_ff -> _:*-_:+) = ^ K , K - cw - ^k* k ac k *k 
v(f s -> k*+k~ ) + r(s° -> _:*-_:+) 1 - AC W - ^k*kC k *k 
r(f s -> __*-__+) - r(_ff -> _:*+_:-) = ^+c^ + ^„Ac^ g 

r(£^ -> K*-K+) + T(B° S ^ K*+K-) \ + AC K * K + szf K * K C K * K 



(5.17) 



(5.18) 

Note that the quantities A k ,± K t here correspond to A s -»+ K -± defined in H37I1 . Therefore, direct CP asymme- 
tries Ax*+k~ an d A K *- K + are determined from the above two equations. Results for various CP -violating 
parameters in the decays 5° — > ^* ± __ :F are shown in Table IXII 
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TABLE XII: CP-averaged branching ratios in B s — ► VV decays (in units of 10 6 ) obtained in various ap- 
proaches. Presented are the pQCD predictions taken from jllll and the QCDF predictions from this work 
and from denoted by BRY. 



Channel 


Class 


QCDF (this work) 


QCDF (BRY) 


pQCD 


Expt [39,64] 


B s - 




T 


—2.8—1.5 


25.2+J-5+J-T 

— — 1 .7 — 3. 1 


—6.2—1.4—1.1 




Bs- 


■+ p°K*° 


c 


. O+2.0+1.7 
—0.6-0.3 


i ^+ 1.0+3.1 
—0.5—1.5 


n +0.09+0. 14+0.00 
—0.07—0.09—0.01 


< 767 


B s - 


-> (oK*° 


c 


i ,+1.5+1.3 
• — 0.5— 0.3 


i ? +0.7+2.3 
—0.3—1.1 


n o, +0.10+0.12+0.07 
-0.07-0.06-0.02 




Bs~ 


-> K*-R*+ 


P 


7 fi +l. 0+2.3 
—1.0—1.8 


g 1+2.5+10.2 
—2.2— 5.9 


7+1.5+3.4+0.5 
—1.2—1.4—0.2 




Bs- 


-> K*°K*° 


P 


6.6+H+H 

" — 1.4— 1.7 


q ,+0.5+11.3 
" — 0.4— 6.8 


7 O+1.9+3.8+0.0 
• — 1.5 — 2.2— 0.0 


< 1681 


B s - 


■+ <pK*° 


p 


rv o 7 +0.06+0.24 
U - J/ -0.05-0.20 


n 4 +0.1+0.5 
u - -0.1-0.3 


n « +0.16+0.27+0. 10 
u ' UJ -0. 13-0. 18-0.04 


< 1013 


B s ~ 


+w 


P 


16 7 +2.6+11.3 
lD./_2 j_ g g 


ni o + l. 1+30.4 
Z1 '°-l. 1-17.0 


-,+8.3+16.7+0.0 
JJ - J -6.9-10.2-0.0 


24.0 ±8.9 


B s - 


■+ 0CO 


P,C 


n iq+0.44+0.47 
u ' io -0.12-0.04 


n 1n +0.05+0.48 
U ' iu -0.03-0.12 


H 1A+0.09+0. 10+0.01 
L '- lo -0.05-0.04- 0.00 




B s - 




Pew 


n i n+0.01+0.09 
u - io -0.01-0.04 


040 +0. 12+0.25 
u,w -0.10-0.04 


n 9 o +0.09+0.03+0.00 
u - ZJ -0.07-0.01-0.01 


< 617 


Bs- 


^p+p- 


ann 


H A8+ 04 + - 73 
u - oo -0.04-0.53 


n o a+0.03+0.60 
u - ,H '-0.03-0.38 


, n +0.2+0.3+0.0 
i - u -0.2-0.2-0.0 




Bs~ 


^p°p° 


ann 


n 0^+0.02+0.36 
u - J ^-0.02-0.26 


n 17 +0.0l+0.30 
u ^ ' -0.01-0.19 


n <ri +0.12+0.17+0.01 
U - Ji -0.11-0.10-0.01 


< 320 


Bs- 


^p°co 


ann 


004 +0 0+0 005 
"■"^-0.0-0.003 


<0.01 


n nm+ 002 +o ooi+o.ooo 

u ' uu/ -0.001-0.001-0.000 




Bs- 


-> coco 


ann 


n iQ+0.02+0.21 
u - 1 ^-0.02-0.15 


n i i +0.01+0.20 
u - 11 -0.01-0.12 


f> ™+0.09+0.13+0.01 
u "^-0.08-0.07-0.00 





VI. B s -» VV DECAYS 
A. Branching fractions 

In two-body decays B u>d —> PP,VP,VV, we have the pattern VV > PV > VP > PP for the branching 
fractions of tree-dominated modes and PP > PV ~ VV > VP for penguin-dominated ones, where the fac- 
torizable amplitude for B — ► VP(PV) here is given by (V(P)|7 jU |B)(P(V)|/' i |0). The first hierarchy is due to 
the difference of decay constants fy > fp and the second hierarchy stems from the fact that the penguin am- 
plitudes are proportional to + r^a^, #4 + r v x a^, a\ — r^ci(, + r^af,, respectively, for B —> PP, PV, VP, VV 
with r£ ~ 3> r%. The same is also true in the B s sector. From Tables ITTTl I VIII and IXIII we find 

m(B s -> p~k*+) > m(B s -» p-K + ) > m(B s -> %-k*+) > -> 7zta: + ), 

SB[B S -> tf+JT) > -> K*~K + ) ~ -> K* + K~) > SB[B S -> £*-.fi:* + ), (6.1) 

for tree- and penguin-dominated 5 V decays, respectively 

There exist two QCDF calculations of B 4 — ► VV Q, @]. However, only the longitudinal polarization 
states of 5 5 — > VV were considered in |0]. The analysis in this work differs from Beneke, Rohrer and Yang 
(BRY) J^l] mainly in three places: (i) the choice of form factors, (ii) the values of the parameters and 
c/»a, and (iii) the treatment of penguin annihilation contributions characterized by the parameters j8,- [see 
Eq. (12.101 )1 for penguin-dominated VV modes. First, the form factors for B s — ► K* and B s — > transitions 
we employ in Eq. (13.61 ) are smaller than the ones (13.101 ) used by BRY. Second, BRY applied the values 
p A (K*(j>) = 0.6 and M K *<P) = - 40 ° obtained from a fit to the data of B -> to study £ -> and 
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B s — ► VV decays. However, as pointed out in [43], the parameters p A (K*p) ps 0.78 and <f> A (K*p) « -43° fit 
to the data of B — ► iPp decays are slightly different from the ones p A {K*ty) and <p A (K*<p) . Therefore, within 
the framework of QCDF, one cannot account for all charmless B — ► W data by a universal set of p A and 0a 
parameters. This explains why the B — ► ^T*p branching fractions obtained by BRY are systematically below 
the measurements. In this work, we choose p A = 0.70 and <p A = —55° (cf. Table [Tib to describe B s — > W 



decays. Third, as noticed in |43|], there are sign errors in the expressions of the annihilation terms A3' and 
A 3 '° obtained by BRY. As a consequence, BRY claimed (wrongly) that the longitudinal penguin annihilation 
amplitude j8 3 ° is strongly suppressed, while the /$ 3 ~ term receives sizable penguin annihilation contribution. 
This will affect the decay rates and longitudinal polarization fractions in some of B — > K*p modes, as 
discussed in details in [43]. In spite of the above-mentioned three major differences in the calculations of 



this work and BRY, it turns out that the calculated rates and fi shown in Tables IXTTl and IXTVl respectively, 
are similar for most of the B s — > VV modes. 

Recently CDF has reported a new measurement of B s 

m(B s 00) 



00 m 

1.78 ±0.14 ±0.20) x 10~ 2 . 



(6.2) 



Using the branching fraction of B s 
to 



J/Y<p from PDG 13911 . updated to current values of f s /fd, this leads 



^(B,^00) = (24.0±2.1±2.7±8.2) x 10 



-6 



(6.3) 



where the error is dominated by the last uncertainty coming from the J/yfcj) branching fraction error. This 
new measurement is slightly larger than the previous one of (14ly) x 10~ 6 16511 . Our prediction £8{B S — > 
00) 16.7 x 10~ 6 is consistent with experiment. 

A few words on the penguin-dominated decays B s — ► <f>K*° and B s — > ft)0. Their branching fractions of 
order 10 -7 are much smaller than other penguin-dominated K*K* and 00 modes. This is because B s — > <j)K* Q 
is induced by the b —> d penguin transition. The amplitude of B s — > ft)0 reads 



V2A, 



d P" a ^ +2(X 3+ 2 <EW 



(6.4) 



The branching fraction due to the QCD penguin (X3 =03 + ^5 is small, only at the level of 10~ 7 . Moreover, 
there is a partial cancellation between QCD and electroweak penguin contributions, making its rate even 
smaller. 6 

As seen from Table IXIII pQCD predictions for the color-suppressed tree-dominated modes p°K*° and 
coK* are much smaller than the QCDF results, whereas SB{B S — > 00) = ^(35 x 10~ 6 ) is much larger than 
QCDF and the CDF measurement 16511 . 

In analog to Eq. (14. 12b . there are three SU(3) relations relating the rates of B s — > VV and Bj — > VV: 



?(B s ^K* + p-)^ 



i(Bd P + P 
3S{BT 



m{B s ^K* + K*-) K,m{B a 



K*-p + ) 



K*°p-) 



(6.5) 



6 It was argued in |)8|] that the color-suppressed tree amplitude <Xi is the largest partial amplitude in the decay By — > 0)0 . 
We found that this decay is still dominated by the QCD penguin, though the contribution from a-i is not negligible. 
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TABLE XIII: Same as Table IXHl except for direct CP asymmetries (in %) in the B s — > VV decays. 



Channel 


Class 


QCDF (this work) 


QCDF (BRY) 


pQCD 


S s - 




T 


1 ,+1+4 

ii -i-i 


T+l+2 

J -l-3 


o 9 + 1.0+1.2+0.4 
°' z -1.2-1.7-l.l 


Bs~ 


+ p K *Q 


C 


17-25 


97+5+34 

z -7-27 


M 0+3.2+17.1+4.4 
oi -°-4.7-22.8-2.3 


B s - 


-> (0K*° 


C 


jrrv+20+21 
JU -15- 6 


cm+10+31 
■>'*- 7-43 


-r 9 1+4.8+19.7+5.5 
OZ.l j 912 61 9 


B s ~ 


-> K* K*+ 


P 


71+1+2 


9+0+40 
z -0-15 


q 7+0.4+3.3+0.3 
*- J -0.7-3.6-0.2 


Bs- 


-> K*°K*° 


P 


n 4+O.8+O.6 
u - H -0.5-0.4 


1 +0+1 
i -0-0 





B s - 




P 


q+3+4 
y -l-6 


i 7+4+9 
1 '-5-9 





B s - 




P 


n 9 +0.4+0.5 
U ' Z -0.3-0.2 


i+0+l 
i -0-0 





Bs- 


-y 00) 


P,C 


O+3+20 
~ 8 -l-15 


8+3+102 
8 -3- 56 


i f. +0.6+2.4+0.6 
J - u -0.6-2.4-0.2 


Bs- 




Few 


87 + 1 + 10 
8:S -0-36 


1Q+5+56 
19 -5-67 


10 1+0.9+1.6+1.3 
iu - i -0.9-l. 8-0.5 


B s - 




ann 







9 1+0.2+1.7+0.1 
^" 1 -0.1-1.3-0.1 


Bs ~ 


+ p°p° 


ann 







9 1+0.2+1.7+0.1 
1 -0.1-1.3-0.1 


B s - 




ann 







n +0.7+2.7+1.0 
°- u -0.5-3.9-0.4 


B s - 


-> coco 


ann 







9 n +0.1+1.7+0.1 
^• u -0.1- 1.3-0.1 



Numerically, we have 

2l.6tH=24.2tli, 7.4^=8.9+^, 6.6 ± 2.2 = 9.2 ± 1.5 (6.6) 
in units of 10~ 6 , where use of the theoretical calculation of ^(fij — > K*~ p + ) from J43I1 has been made. 



B. Direct CP violation 



Direct CP asymmetries in QCDF and pQCD approaches are summarized in Table IXIIll 



C. Polarization fractions 



For charmless B — > VV decays, it is naively expected that the helicity amplitudes ^4 (helicities h = 
0, — , + ) for both tree- and penguin-dominated B — > VV respect the hierarchy pattern 

^:^:^ = l:(^y.f— Y- (6.7) 
\ m b J \ m b ) 

Hence, they are dominated by the longitudinal polarization states and satisfy the scaling law, namely [66], 

/^w,=^4y £=i+*(-v (6.8) 

with ft,f±, f\\ and / T being the longitudinal, perpendicular, parallel and transverse polarization fractions, 
respectively, defined as 

f = Fa \^a\ 2 (( _ qx 

U - T " |^ |2 + |^|2 + |^ ± |2' 
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with a = L, ||, X. In sharp contrast to the pp case, the large fraction of transverse polarization of order 
0.5 observed in B — > K*p and B — > K*ty decays at B factories is thus a surprise and poses an interesting 
challenge for any theoretical interpretation. Therefore, in order to obtain a large transverse polarization in 
B —> K*p,K*<p, this scaling law must be circumvented in one way or another. 

As pointed out by Yang and one of us (HYC) |43|], in the presence of NLO nonfactorizable corrections 
e.g. vertex, penguin and hard spectator scattering contributions, effective Wilson coefficients d\ are helicity 
dependent. Although the factorizable helicity amplitudes X°, X and X + defined by Eq. (12.41 ) respect the 
scaling law (16.71 ) with Aqcd/W replaced by 2my jm^ for the light vector meson production, one needs to 
consider the effects of helicity-dependent Wilson coefficients: =sz/~ = f(aj)X~ /[f(af)X ]. For some 
penguin-dominated modes, the constructive (destructive) interference in the negative-helicity (longitudinal- 
helicity) amplitude of the B — ► VV decay will render f{a[) S> f{a®) so that si~ is comparable to g/° 
and the transverse polarization is enhanced. For example, fa{K p ) ~ 0.91 is predicted in the absence of 
NLO corrections. When NLO effects are turned on, their corrections on aj will render the negative helicity 
amplitude stf~ (B° —* K*°p°) comparable to the longitudinal one £/°(B° — > K*°p°) so that even at the short- 
distance level, fa for B° — ► K*°p° can be as low as 50%. However, this does not mean that the polarization 
anomaly is resolved. This is because the calculations based on naive factorization often predict too small 
rates for penguin-dominated B — ► VV decays, e.g. B — > K*ty and B — ► K*p, by a factor of 2 ~ 3. Obviously, 
it does not make sense to compare theory with experiment for faj as the definition of polarization fractions 
depends on the partial rate and hence the prediction can be easily off by a factor of 2 ~ 3. Thus, the first 
important task is to have some mechanism to bring up the rates. While the QCD factorization approach 
relies on penguin annihilation [66], soft-collinear effective theory invokes charming pengu in Il67ll and the 



final-state interaction model considers final-state rescattering of intermediate charm states B681. 1691 . 17011 . A 



69, 



nice feature of the (S — P)(S + P) penguin annihilation is that it contributes to £/ Q and with similar 
amount. This together with the NLO corrections will lead to fi ~ 0.5 for penguin-dominated VV modes. 
Hence, within the framework of QCDF we shall assume weak annihilation to account for the discrepancy 
between theory and experiment, and fit the existing data of branching fractions and fi simultaneously by 
adjusting the parameters p A and <f> A . Then using this set of annihilation parameters as a guideline, we can 
proceed to predict the rates and fa for other VV decays of the B u j s mesons. 

The longitudinal polarization fractions in B s — > VV decays obtained in the QCDF and pQCD approaches 
are summarized in Table IXIVI Transverse polarization effects are sizable in penguin-dominated B s — * VV 
as expected. However, the pQCD calculations indicate that fa ~ fa ~ I even for the color-suppressed tree- 
dominated decays B s — > K*°(p ,co). This is an astonishing result and should be checked by experiment. 
Polarization fractions of B s —> <jxp will be studied soon by CDF. It will be very interesting to see if the 
transverse polarization is also important in the penguin dominated B s decays. 
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TABLE XIV: Same as Table IXIll except for the longitudinal polarization fractions in the B s — > V V decays. 



Channel 



Class 



QCDF (this work) 



QCDF (BRY) 



pQCD 



B s 
Bs 
B s 
Bs 
B s 
B s 
B s 
B s 
Bs 
Bs 
Bs 
Bs 
B s 



p-K*+ 
. p0jr o 

coK*° 

■ K*K*^ 

K*°r° 

(j>K* 
<pco 

■<pp° 

■ p+p- 
p°p° 

P °(0 

coco 



T 
C 
C 
P 
P 
P 
P 

PC 

Pew 
ann 
ann 
ann 
ann 



0.92 
0.90; 

0.90; 

0.52; 
0.56; 
0.43; 
0.36 
0.95 
0.88 



-0.01+0.01 
-0.02-0.03 
-0.04+0.03 
-0.05-0.23 
-0.03+0.03 
-0.04-0.23 
-0.03+0.20 
-0.05-0.21 
-0.04+0.22 
-0.07-0.26 
-0.02+0.21 
-0.02-0.18 
-0.03+0.23 
-0.04-0.18 
-0.01+0.00 
-0.02-0.42 

-0.01+0.02 
-0.00-0.18 

1 
1 
1 
1 



0.921°, 
0.93+° 
0.93+°, 
0.67+° 
0.63+° 
0.40+°,; 
0-43+°; 



01+0.05 
01-0.08 
02+0.05 

03- 0.54 
02+0.05 

04- 0.49 
04+0.31 

05- 0.26 
00+0.42 

00- 0.29 
01+0.67 

01- 0.35 

00+0.01 
00-0.34 



rv oi +0.03+0.09 
u - 01 -0.04-0.12 



0.937 



-0 



0.4551°, 
0.532+° 



0.438+° 
0.497+° 
0.712+° 



0.619 



0.4431° 
0.8701°, 



001+0 

002- 
004+0 

003- 
003+0 
002-0 
051+0 
040-0 
057+0 
048-0 
032+0 
030-0 
036+0 
032-0 
000+0 
075-0 
002+0 
002-0 



002+0.000 
003-0.002 
069+0.006 
.043-0.009 
.035+0.023 
.029-0.013 
.021+0.037 
.023-0.015 
.006+0.000 
.038-0.000 
.027+0.000 
.037-0.000 
.025+0.000 
033-0.000 
054+0.009 
061-0.004 

.009+0.009 
.003-0.004 



TABLE XV: Direct CP asymmetries (in %) in B s — ► VV decays via £/-spin symmetry. 



Modes 



?(10" 6 ) 



A CP (%) 


Modes 


A C p(%)(U-spm) 


A CP (%)(QCDF) 


09 + 1+ 5 
JZ -3-24 


B°s~ 


->K*+p- 


-10.2 


-lilt 


i r+4+16 
iJ -8-14 


t~ 


■+K*°p° 


42.3 


46+18 


_ 4 +o+3 

+-0-3 


t- 


^K*+K*- 


18.7 


21+1 


_ 14 +l+6 




t- 
t~ 


^p+p- 


0.5 



4+1.0 

U -^-0.6 





t 



K*-p+ 

K*°p° 

p+p- 

K*+K*- 



1.1+4.8 
-1.0-5.5 
0.6+3.5 
-0.5-3.5 
+ 1.5+2.4 
-2.6-1.5 
fi +0. 1+0.2 
u -°-0.1-0.3 



8.9 
4.6 
25.5 



0.15 



+0.02+0.11 
-0.01-0.12 



D. U -spin symmetry 

Analogous to the B s — > PP sector, U -spin symmetry leads to the following relations: 

A CP (B s ^K* + p-) = -A CP (B d ^K*-p + ) 

A CP (B S ^K* + K*-) = -A CP (B d ^p + p 



*+„^ . „ ts ^*-„ + ,^(Bd^K*-p+) x{B s ) 



(B s ^K*+p-) z(B d y 
Sg{B d ^p+p-) t(B s ) 



A CP (B S ^K*°K*°) = -A CP (B d -> K*°K*°) fffi - ^"f °j 4^4, (6-10) 

a f^rW- a (P. : w ^^^V) < B °) 

A CP (B S ^K p ) - -A CP (B d ^K p ) ^ - ^ 0p0) ^ 

AcpiB^p+p-) = -A Cp{ B d ^K n ^ m%m 



(B s ^p+p-) z(B d y 

In Table lXVl we compare the results of CP asymmetries inferred from U -spin relations with the direct QCDF 
calculations. It appears that U -spin symmetry works well in the VV sector. 



32 



Assuming that the transverse amplitude can be expressed as a single dominant contribution which may 
arise from new physics, £/-spin symmetry implies that the transverse amplitudes of B s — > VV can be related 
to the U -spin related decays in the B d sector via 117 ill 



f T (B d ^K*°K*°) 



Vts 



v td 



hi 



Therefore, 



f T {B s -> K*°K*°) 
f T (B d - K*°K*°) 
f T (B d -> $K*) 
f T (B s ^^K*°) 



(25.5 ±6.5) 
(19.3±4.9) 



^ T (B d ^$K*°) 
£/ T (B d ^$K*°) 

<%{B d 



Vu 



Vtd 



hi 
h s 



(6.11) 



m{B, 

3&{B S 



(j>K*°) 



®(B d 

The polarization measurement in the B d decay thus allows one to predict the transverse polarization in the 
B s decay. 7 Using the data H48TI 

x 1(T 6 

-6 



(1.28+^) x 10 



®(B d 

38{B d (j)K*°) = (9.8 ± 0.7) x 10 
and QCDF predictions for 3{B S -> K*°K*°) and 3B(B S 
f T {B s ^ K*°K*°) = 1.02 ±0.28, 



<j)K* Q ) 



(6.12) 



o.so±8:ii, 



h(Bd 

f L (B d ^ (j>K*°) = 0.48 ±0.03, 
<j)K*°), we obtain 



f T {B s ^(j)K* ) = 0.73 ±0.19. 



(6.13) 



(6.14) 



It is obvious that the central value of the predicted fr(B s — ► ^*°^*°) via £/-spin symmetry is too large. 
Note that there is a discrepancy between the QCDF prediction of k*°K*°) = (0.6±g;|) x 10~ 6 |43|] 

Jd ^ J £*°^*°) = (1.28+°;^)-' 



and the BaBar measurement SB{B d -> /T ^* ) = (1.28+HZ) x 10~ 6 ILZJ. We need to await a more precise 



measurement of B d — ► k*°K*° in order to have a more accurate prediction of its transverse polarization 
fraction via t/ -spin symmetry. 



E. Time-dependent CP violation 

In principle, one can study time-dependent CP asymmetries for each helicity component, 

r(fl{t)^v h vn+rw(t)-+v h vj i ) 

= Sh sin (Am s t) — Ch cos (Am/), (6.15) 

where the effects of the width difference of the fi v mesons have been neglected. From Table IXIII we see 
that there is only one decay mode of particular interest, namely, B s — ► <f)<f). Indeed, this could be the most 
promising channel for the forthcoming LHCb experiment. This channel is a pure b — > sss penguin-induced 
process and hence provides an ideal place for exploring the signal of New Physics via B s — B s mixing and/or 
the penguin process. The other decays such as B s — > pp,p°co, (0(0 proceed through weak annihilation. The 



7 Based on SU(3) flavor symmetry, it has been shown in M72I1 that the transverse polarizations of B s — > 00 and 



B s -> 0/T° can be related to B d -> 0^*° and B d -> K*°K*°, respectively. 
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modes 0ft) and 0p receive QCD penguin and electroweak penguin contributions, respectively, but their 
rates are too small. A straightforward calculation gives 

^l = (5.9+^)x10- 6 , C L = (-0.5t Q Q t\ A 5 )%, S L = (-05t° Q t\i)%, (6.16) 

for the longitudinal component of B s —> 00. Note that Sl is found to be positive and small < 0.02 in |@], 
while our result is negative for Sl- An observation of large CP violation in this decay will rule out the 
scenario of minimal flavor violation. Time-dependent CP violation will be studied at LHC. If LHCb is 
upgraded to accumulate data sample of 100fb _1 , the sensitivity of Sb s ^^ l will reach the level of 0.01 ~ 
0.02. 



VII. CONCLUSIONS 

We have re-examined the branching fractions and CP-violating asymmetries of charmless B s — > 
PP, VP, VV decays in the framework of QCD factorization. We have included subleading power corrections 
to the penguin annihilation topology and to color-suppressed tree amplitudes that are crucial for resolving 
the CP puzzles and rate deficit problems with penguin-dominated two-body decays and color-suppressed 
tree-dominated 71° 71° and p°7r° modes in the B u ^ sector. Our main results are: 

i) . Many model-independent relations for CP asymmetries and branching fractions of B c \ and B s decays 

can be derived under U -spin and SU(3) symmetries for PP, VP, VV modes. In general, they are either 
experimentally verified or theoretically satisfied. There are also a few U -spin relations for transverse 
polarizations in B s — > VV decays. 

ii) . For the B s — > K transition form factor, we use a smaller one, F B " K 0.24 at q 1 = obtained by the 

lattice calculation, to avoid too large rates for B s — > K + n ,K + K decays. 

iii) . Both QCDF and SCET indicate that the penguin-dominated decay B s — > r\'r\' , the analog of B —> Kr\' 

in the B s sector, has the largest branching fraction of order ~ 50 x 10~ 6 in two-body hadronic decays 
of the B s meson, whereas the pQCD approach claims that S3(B S — > ITH') ~ 35 x 10~ 6 is the largest 
one. 

iv) . Even at the decay rate level, there are some noticeable differences between various approaches. The 

branching fractions of the color-suppressed tree-dominated decays obtained by pQCD, for exam- 
ple, B s -> K n ,K°ri^,K* Q n ,p K ,o)K ,K* ri' are typically smaller by one order of magnitude 
than that of QCDF and SCET. For example, 3${B S -> p°K°) is predicted to be of order 1.9 x 10~ 6 
by QCDF, but it is only about 0.08 x 10~ 6 in pQCD. In the QCDF approach, many of the above- 
mentioned decays get a substantial enhancement from the power corrections to the color-suppressed 
tree topology. 

v) . The decay rate of B s — > §r\' is sensitive to the 5 S — > (j) transition form factor A Q (0). For A Q (0) = 

0.474 obtained by the light-cone sum rule method, a near cancelation between B s — > 0T] S . and B s — > 
(j>T] q occurs in the decays B s — > 0T]', so that its branching fraction, of order 10~ 7 , becomes very 
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small. However, if the value A ( , (0) = 0.30 favored by many other model calculations is employed, 
then B s — ► 07],, and B s — > 07]^ will contribute constructively to 5 A — > 07]' so that ^(B s — > 07]') = 
2.2 x 10~ 6 and £§{B S — > 07]) = 1.0 X 10~ 6 . Hence, it is very important to measure the branching 
fractions of B s — ► 07] to gain the information on the form factor Aq 1 ^. 

vi) . Measurements of CP-violating asymmetries can be used to discriminate between QCDF, pQCD and 

SCET approaches: 

(a) Both QCDF and pQCD predict a positive sign for A CP (B S -> K°7t°), whereas SCET leads to a 
negative one. This can be traced back to fact that Acp(Bd —> K°7i°) is positive in SCET, while it 
is negative inferred from the CP-asymmetry sum rule, SU(3) relation and the topological quark 
diagram analysis. 

(b) For color-suppressed tree-dominated decays B s — ► K* Q n ,p°K , coK ,K*°ri', QCDF and pQCD 
results are of the same sign, whereas SCET predicts opposite signs for these modes. In the 
QCDF approach, the signs of these CP asymmetries are governed by the soft corrections to 
Since the corresponding rates of these decays are very small in pQCD, as a consequence, the 
CP- violating asymmetries predicted by pQCD are very large, of order 0.50 or even bigger. 

(c) In the QCDF framework, the penguin-dominated decays B s ->■ K°(j),K* K ,K*°K° have non- 
vanishing CP asymmetries, though very small for the last two modes, whereas leading order 
pQCD predicts no CP violation for these three decays. 

vii) . Mixing-induced CP asymmetries of the penguin-dominated decays B s — > K°K°, 7] ^7]^, 07]', 00 are 

predicted to be very small in the SM. Especially, we found Sg s _,^ L ~ —0.5%. They are sensitive to 
New Physics and provide possibilities of new discoveries. While both QCDF and pQCD approaches 
predict S Ss ^ Ks<tl ~ ^(0.70), the SCET result of 0.09 or -0.13 is dramatically different. 

viii) . Due to soft power corrections to the color-suppressed tree amplitude, we find that such effects will 

convert the sign of mixing-induced CP violation 5/ into the positive one for the color-suppressed 
decays B s — ► Ks(n { \ 7], 7]'). Therefore, even the measurements of the sign of Sg s _^. Ks t^o ^ ^ will be 
helpful to test if "a2" has a large magnitude and strong phase. 

ix). Transverse polarization effects are sizable in penguin-dominated B s — > VV as expected. However, 
the pQCD approach predicts that fa ~ fj ~ A even for the color-suppressed tree-dominated decays 
B s — > K*°(p°, ft)). This should be tested by experiment. 
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